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THEORETICAL STUDY OF THE USE OF VARIABLE GEOMETRY 


IN THE DESIGN OF MINIMAL -CORRECTION 
V/STOL WIND TUNNELS 

By Harry H. Hey son 
Langley Research Center 

SUMMARY 

Theoretical study indicates that, if either the width -height ratio or the model height 
is properly scheduled, the vertical interference velocities at the lifting system can be 
reduced to zero in closed -on -bottom -only tunnels. Reductions in interference at the tail, 
nonuniformity of interferences, and minimum speed for recirculation -free testing can be 
obtained simultaneously; however, these reductions are much greater in the case of vari- 
able width-height ratio. Variable width-height-ratio operation of a closed tunnel can 
reduce the interference at the lifting system by a factor of 2 or 3. This configuration, at 
low speeds, can reduce interference at the tail, nonuniformity of interference, and mini- 
mum speed for recirculation-free testing to almost negligible values. 

INTRODUCTION 

The large magnitude of wind-tunnel wall interference in V/STOL testing has been 
documented for several years (for example, ref. 1). The magnitude of these wall effects 
is calculable by theory (refs. 2 to 4), and substantial experimental evidence indicates that 
the theoretical results can be used to correct wind-tunnel data for wall effects (refs. 5 
to 11). On the other hand, there is little comfort to be found in applying theoretical cor- 
rections for large wall-induced angularities (which may exceed 10°) when it is recognized 
that large corrections, with the associated large flow distortions, may effectively alter 
the configuration of the model from a viewpoint of aerodynamic equivalence. Further- 
more, omission in the corrections, whether by simplification or oversight, of factors 
which influence the corrections (for example, finite span, load distribution, changes in 
model or tail position as a function of angle of attack (ref. 12) may produce significant 
errors when the corrections are large. Thus, it is highly desirable to find some means 
whereby the magnitude of V/STOL wall effects can be significantly reduced; however, 
unless the wall effects are actually reduced to zero, it is still required that a means be 
available to calculate and correct data for such wall effects that remain. 



If all other considerations are equal, wall effects depend directly upon the ratio 
between an area associated with the model and the cross-sectional area of the wind-tunnel 
test section. Thus, the classic solution to the problem of large wall effects has been to 
either rim the tests in a larger wind tunnel or to reduce the size of the model. While 
this procedure reduces wall interferences, a larger wind tunnel may not be available, 
and there are limits (occasioned by both Reynolds number and the minimum physical 
size of motors and other model components) to the minimum-size model that can be uti- 
lized. In the case of full-scale wind-tunnel tests of an actual aircraft, little or nothing 
can be done to reduce the corrections by means of reducing the area ratio, since the 
model size is fixed and only one or two wind tunnels in the world possess the physical 
size necessary to mount the model. 

Wall-interference theory, in general, obtains only the wall-induced distortions of 
the flow field near the model. In order to correct wind-tunnel data, the effect of the field 
distortion on the model performance must be obtained from the available theoretical 
treatments of the model aerodynamics in free air. The rigorous conversion from inter- 
ference field to correction may prove difficult even in the case of simple conventional 
airplanes for which 50 years of theoretical studies provide a vast body of material with 
which to work. For many V/STOL models, essentially no theoretical background exists. 
Under such circumstances, it may prove impossible to use the wall-interference calcula- 
tions to provide rigorous corrections. 

The alternative is to develop test-section configurations with zero or near-zero 
wall effects. Several such configurations have been found and utilized in conventional 
testing where the wake essentially passes directly downstream (refs. 13 and 14). In addi- 
tion, studies (such as ref. 15) of slotted wind tunnels developed primarily from the con- 
sideration of transonic choking, have disclosed numerous other configurations which 
possess zero, or near zero, "lift interference" for the case of an undeflected wake. At 
the present time, a number of slotted tunnels intended for low -speed V/STOL testing are 
under construction. Wall interference in these tunnels is under study, both theoretically 
and experimentally (refs. 16 to 20). Theoretically, it should be noted that zero longitud- 
inal interference is not obtained simultaneously with zero vertical interference at a single 
slot opening under undeflected -wake conditions (ref. 15). It is not reasonable to assume 
that this trend will be changed by the inclusion of large wake deflections; indeed, experi- 
mental observations (refs. 10, 19, and 20) indicate that the two interferences do not vanish 
simultaneously in V/STOL testing. Slots markedly reduce the wall interference in all 
cases tested. The reduction is not completely to zero, however, and the current lack of 
an adequate theory to account for the remaining interference when testing with highly 
deflected wakes poses a certain element of risk in the use of slotted tunnels. 
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References 10, 19, and 20 report tests in which a closed tunnel is completely 
opened behind a given longitudinal station in the test section. This truncated test section 
yielded very small corrections when the opening commenced at the model location. Such 
results might be anticipated from the theoretical study presented in reference 21. It will 
be noted, however, that references 19 and 21 indicate large longitudinal gradients of inter- 
ference. Thus, the truncated test section would be expected to show large wall interfer- 
ences on pitching moments and even on the performance of very long lifting systems 
(such as tandem -rotor helicopters). In addition, the theoretical treatment of wall inter- 
ference in such sections is exceedingly difficult (ref. 21) even for an undeflected wake. 

The streamline matching technique of reference 22 appears feasible and is pres- 
ently under study. In this technique, the free-air flow field of the model is calculated at 
the location of the walls, and the flow angles at the walls are then adjusted (by means of 
louvers, air jets, etc.) to coincide with the calculated free-air field. Theoretically, this 
technique totally removes the wall interference. The calculation of the flow field for an 
arbitrary high-lift model is difficult and may require several hours on a large digital 
computer. It will be noted that the calculation is, in general, a far -field calculation 
(allowing some simplification) except at the regions where the free-air wake is at the 
wall locations. The detail to which the flow in this latter area must be computed in order 
to obtain acceptable results is unknown at present. Experimental studies are underway 
(ref. 22) in order to determine the degree of accuracy in streamline matching that is 
actually required in practice. 

References 2 to 4, which are theoretical treatments of wall effects in more conven- 
tional test sections, note that no single configuration was found which yielded a zero- 
correction test section for all wake -deflection angles (or skew angles). On the other 
hand, several configurations, such as the closed-on-bottom-only tunnels, were found in 
which the vertical interference due to lift (corresponding to the classical lift interference) 
becomes zero at one wake deflection angle. It would appear that a judicious application 
of variable -geometry to such configurations, with the geometry alterations applied as a 
function of wake deflection, would result in a wind tunnel with zero vertical interference 
due to lift. This single component of interference, in general, produces most of the 
observed wall interference. There would remain much smaller contributions from the 
horizontal interference due to lift as well as from both horizontal and vertical interfer- 
ence due to induced drag. These remaining contributions produce only small effects on 
the model characteristics, and the theoretical corrections can be applied with confidence. 

Reference 9 points out that the deflection of the rolled-up wake should be used in 
applying the theory of references 2 to 4. This deflection is in contrast to the original 
usage of the deflection as obtained directly from momentum theory. Since the rolled-up 
wake has only approximately half the momentum deflection (at the center of lift), it is 
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necessary to provide only sufficient variable geometry to null the correction over a much 
smaller range of wake deflections. This feature alone significantly reduces the severity 
of the required changes. 

Several classes of wind tunnels are considered herein. These include closed-on- 
bottom -only configurations in which either the model height or the tunnel width -height 
ratio is varied as a function of wake deflection. Similar closed-tunnel arrangements 
are also considered since, even though the interference factor cannot be nulled in these 
tunnels, it can be minimized. 

Since the corrections are conditioned by the model size and configuration, calcu- 
lated results are presented for straight wings of various span-width ratios as well as 
for a vanishingly small lifting system. The effect of configuration is indicated by calcu- 
lated results for wings of differing sweep and for a lifting rotor of finite size. 

The corrections at the model tail are also considered since these corrections differ 
from those at the lifting system. Calculations for several different tail lengths are made 
and presented. The variation of interference with tail length provides an indication of 
the magnitude of the interference gradients at the lifting system itself. These latter gra- 
dients are significant since they can produce significant effects on wind-tunnel data, par- 
ticularly with regard to pitching moment (ref. 18). 

The lateral distribution of interference over wings with finite span-width ratio in 
the variable -geometry tunnels is studied briefly. Comparisons are made between the 
distributions in the fixed and variable-geometry tunnels. 

One of the most significant developments in V/STOL testing has been the recent 
discovery by Rae of the University of Washington (ref. 23, see also refs. 18 and 24) of the 
effect of recirculation on limiting the minimum speed at which acceptable data can be 
obtained in a closed wind tunnel. Similar effects may be expected in the closed-on- 
bottom -only configuration as well, if for no other reason than the presence of a boundary 
layer on the closed floor of the wind tunnel (ref. 25). In any event, the occurrence of 
this phenomenon can be related to the downstream distance of the intersection between 
the theoretically straight wake and the floor. Thus, this effect is studied herein by 
examining the relative behavior of the foregoing impingement distance in the various 
wind-tunnel configurations. 

The provision of variable geometry in a wind tunnel poses many related practical 
problems in the design of the entrance, diffuser, balance, fan, and drive systems. Sev- 
eral of these problems are discussed briefly in connection with each of the tunnel con- 
figurations examined herein. In certain types of variable -geometry tunnel, the same fea- 
tures used to negate V/STOL wall interference can be used to provide test capabilities 
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other than those usually available in a fixed geometry wind tunnel. Several of these addi- 
tional capabilities are noted. 


SYMBOLS 

The selection of a single set of symbols and definitions for the wide variety of aero- 
dynamic systems encountered in wind-tunnel work does not allow complete conformity 
with existing practice. The following list sets forth the terminology used herein. Posi- 
tive directions are self-consistent, that is, all forces, directions, and velocities are posi- 
tive when directed in the positive sense of the chosen axes, as illustrated in the sketch 
below. Similarly, all moments and angles are positive in the direction of the right-hand 
rule with the chosen axes. Certain unusual features result (e.g., a negative induced 
velocity Wq results from a positive lift L). The reader should carefully consider the 
following definitions and make appropriate conversions for his own use. It should be 
noted that these definitions as applied to wings differ from reference 12 in that the origin 
is chosen to be at the aerodynamic center in the present paper. 



A m momentum area of lifting system 

cross-sectional area of test section 

b distance from right-hand side of test section (viewed from behind) to origin 

in model 
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B 

C L 

COB 

Di 

h 

H 

h 

L 

R 

s 

u 

u 0 

V 

w 

w 0 

w h 

x f 


semiwidth of test section 
lift coefficient 

abbreviation for closed on bottom only 
induced drag 

height of model aerodynamic center above test- section floor 
semiheight of test section 

tail length, distance of tail aerodynamic center behind model aerodynamic 
center at zero angle of attack 

lift 

rotor radius 
semispan of wing 

interference velocity in x-direction 

mean, or momentum theory, value of model induced velocity along X-axis, 
positive rearward 

tunnel airstream velocity 

interference velocity in z -direction 

mean, or momentum theory, value of model induced velocity along Z-axis, 
positive upward 

hovering induced velocity, value of Wq when hovering with both zero 
forward speed and drag 

theoretical impingement distance of straight-line wake on floor measured 
from model 
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x,y,z 


value of xj when model is fixed at center of the tunnel 

distances from origin along X,Y,Z axes, positive when in positive direction 
of axis 


X,Y,Z Cartesian axes with origin at the model aerodynamic center, parallel to tunnel 
axes; X positive rearward, Y positive to right when viewed from behind, 
Z positive upward 


a 


angle of attack 


u,D 


width-height ratio of tunnel, B/H 


interference factor for horizontal interference velocity due to drag, defined 

A m 

implicitly by Au D = 6 n n - — u 


%D u o 



interference factor for horizontal interference velocity due to lift, defined 

A m 

implicitly by Au L = 6 u L — w Q 


6 w,D 


interference factor for vertical interference velocity due to drag, defined 

A m 

implicitly by Aw D = 6 W D — u Q 


w,L 


interference factor for vertical interference velocity due to lift, defined 
implicitly by Aw L = 6 WjL ^ w 0 


AH height of model aerodynamic center or rotor center above wind-tunnel 

centerline 


Au D horizontal interference velocity due to drag 

Au L horizontal interference velocity due to lift 

Aw total vertical interference velocity 

Awp vertical interference velocity due to drag 
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vertical interference velocity due to lift 


Aw l 

AS w l difference between interference factors for vertical interference velocity due 
to lift at the tail and at the center of lift 


£ ratio of test-section semiheight to height of model aerodynamic center above 

floor, H/h 

A wing sweep angle, angle between lateral axis of model and lifting line, 

positive rearward 

a span -width or diameter -width ratio, s/B or R/B 

a. span-width ratio of tail, ™ semispan 

t B 

X effective wake skew angle, angle between effective center of rolled-up wake 

and negative Z-axis, positive rearward 

X m momentum -theory value of wake skew angle at the lifting system 


Superscript: 


referenced to high-speed section with y = 2.0 (for example 6 T = 6 T — - 

W,1j W,li a t 

* 

so that Aw l = 6 w?l — w ( 


RESULTS AND DISCUSSION 


General Observations 


Closed tunnels .- Figure 1 displays the interference factors in closed tunnels for the 
vertical interference velocity due to lift 6 W ^ calculated for a zero-span model (refs. 4 
and 12) as a function of the effective wake skew angle x- Figure 1(a) compares the inter- 
ference factors for tunnels of different width-height ratio y. Figure 1(b) presents a 
comparison of the interference factors for different heights above the centerline AH in 
a tunnel with y = 1.5. The results presented in figure 1(b) were obtained from the com- 
puter programs described in reference 12 by noting that 


AH _ 1 « 

H C 


( 1 ) 
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Irrespective of the combination of vertical model position AH/H or tunnel width - 
height ratio y chosen, the interference factor 6 W ^ is always negative. This fact 
indicates the presence of an upwash at the model, which is expected since the completely 
closed tunnel, by restricting the normal downwash required for lift, effectively presents 
an upwash to the model. It is obvious that no scheme of altering the wind-tunnel propor- 
tions or model position can result in a completely closed tunnel with zero interference. 
Nonetheless, it may be possible to reduce the interference somewhat by a judicious choice 
of these parameters as a function of wake angle. Since, at the outset, the degree of 
reduction of interference appears small, discussion of such tunnels will be deferred to 
subsequent sections of this paper. 

Closed-on-bottom-only tunnels .- A similar set of comparisons, again for a zero- 
span model, for the closed-on -bottom-only tunnels is presented in figure 2. In this class 
of tunnels, it is possible to balance the upwash engendered by the restriction of the solid 
floor against the effective downwash produced by the three open boundaries in such a 
manner as to produce either upwash ^6 w ^ < 0) or downwash ^6 W ^ > 0^. In particular, 
it is noted that certain combinations of wake skew angle y either width-height 
ratio y or vertical model position AH/H result in zero vertical interference due to 
lift (o w ^ = 0). This result is not totally unexpected since reference 13 (confirmed by 
ref. 4) has already pointed out that for a small model with an undeflected wake (x = 90°), 
the closed-on-bottom -only tunnel with y - 2 is a "zero -correction" tunnel. Figure 2 
shows clearly that a similar result may be obtained for the deflected -wake case as well, 
provided that either y or AH/H is chosen appropriately for the actual effective wake 
deflection. 

More specifically, if a small computer is tied into the balance system to compute 
the effective wake angle, and then used to adjust y or AH/H to a predetermined 
schedule of values (as a function of y), it should be entirely possible to obtain essentially 
zero boundary interference even for low-speed (high-wake-deflection) tests of V/STOL 
models. Such data would not be totally free of interference effects since there would 
still be horizontal interference velocities and smaller interferences due to drag; however, 
elimination of the vertical interference due to lift would reduce the overall corrections 
to reasonable and acceptable values. 

Closed -On-Bottom -Only Tunnel With Variable Width-Height Ratio 

General concept of configuration .- Schematically, the test section of a closed -on- 
bottom -only tunnel of variable width -height ratio might appear as in figure 3, in which a 
three-dimensional inlet of moderate contraction ratio, say 3 or 4, is followed by a variable 
two-dimensional contraction. The ground plane and the upper and lower lips of the exit 
cone would move with the corresponding members of the entrance cone. 
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The cross-sectional area of the tunnel varies with width-height ratio when this 
ratio is altered by moving only one pair of boundaries. If the upper and lower boundaries 
are moved, the area of the tunnel is greatest at the lowest width-height ratios; if the side 
boundaries are moved, the area is greatest at the largest width-height ratios. Figure 2(a) 
indicates that the lowest speeds (lowest values of y) require low width -height ratios in 
order to negate 6 W j^. As will become evident in subsequent sections, neither the resid- 
ual interferences ^> u j^, 6 W p, and 6 U q) nor the interferences at the tail become zero 
simultaneously with the average value of 6 w L over the primary lifting system. For 
constant lift, the largest values of these interferences should be expected at the lowest 
speeds since Uq, Wq, and are large. It is desirable to have the largest tunnel 

area occur at low speeds, thus accruing the advantages of a low ratio of model area to 
tunnel area for the most severe conditions. Furthermore, at the lower speeds, the 
installed wind-tunnel power is capable of driving a much larger test section than the test 
section provided for at the maximum wind-tunnel speed. Consequently, an increase in 
tunnel cross-sectional area as the speed is reduced tends to provide more efficient use 
of the installed power. Therefore, the arrangement shown in figure 3, where the floor 
and upper boundary move, is preferable to the alternate arrangement where the side 
boundaries are movable. 

The usual form of wind-tunnel speed control is by means of rotational speed changes 
of a fixed pitch fan. This type of speed control is generally relatively sluggish because of 
the large inertias of the heavy rotating equipment and may not be able to cope with the 
more rapid speed changes resulting from the area changes in the tunnel test section. 

Since test programs generally require data at constant speeds, satisfactory operation of 
such a tunnel may require the more rapid action provided by a variable -pitch fan together 
with some form of automatic velocity compensation. 

It should be noted that the additional geometric changes to reduce the tunnel area 
still further and to fold the floor around the stream to form a slotted transonic test sec- 
tion are not so severe as to preclude possible consideration. This additional reduction 
in tunnel area would result in a tunnel of extraordinarily large contraction ratio (on the 
order of 50) which, by itself, should be of considerable advantage in obtaining excellent 
flow quality (low turbulence level, uniform velocity field, and absence of wake shadow). 
Further discussion of this possibility is beyond the scope of this report; however, fig- 
ure 4 illustrates schematically the physical changes that would be required. These 
changes include further reducing the tunnel height, folding the inner surfaces of the side- 
walls inward, folding the ground plane around the smaller stream, and finally opening 
slots in the new closed section. 

Computer programs .- The required schedule of width -height ratio in order to obtain 
§w l ~ 0 (averaged over the lifting system) has been calculated by using the computer 
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programs developed in reference 12. Only minor alterations were required in order to 
have the computer search for the particular value of width-height ratio required for a 
given wake angle and a given wing or rotor configuration. The interference effects at the 
tail are computed simultaneously in order to economize on computer time. In order to 
provide a proper comparison of the effective values of interference, the interference fac- 
tors are altered so as to be referenced at all times to the cross-sectional area of the 
high-speed (y = 2.0) section. Other quantities such as semiheight of the tunnel are also 
referenced to the high-speed-section values. A listing of the complete computer program 
for swept wings is provided in appendix A. A similar listing for rotors is provided in 
appendix B. 


The programs given in the appendixes offer a choice of either of two loadings: for 
wings, the span loading may be either uniform or elliptic; for rotors, the disk-load dis- 
tribution may be either uniform or triangular. All numerical results presented herein 
are for uniform loading, and the interference values presented represent an average over 
the lifting system. Note that results for the variable width-height ratio tunnel are pre- 

A * 

T 

sBiueu in terms uj. o,„ t = o„ T T — in order to include the effect of the changing cross- 

W>1j w,ij Arp ' 

sectional area of the tunnel. 


Required schedule of tunnel width -height ratio for wings .- Figure 5(a) shows the 
required schedule of y to obtain 6 W ^ = 0 for unswept wings of varying span-width 
ratio. The wings are all assumed to be mounted in the center of the tunnel. It may be 
seen that larger values of width-height ratio are required at the higher skew angles as 
the relative span of the model increases. Indeed, an extreme value of y - 3.6 is 
required for a wing with a = 0.75 when the effective skew angle is on the order of 70°. 

Very flat wind tunnels (in excess of y = 2) are undesirable for several reasons. 

In general, the spanwise distribution of boundary interference tends to become substan- 
tially more nonuniform. More significantly, inadequate clearance exists between the tail 
and the floor at large angles of attack, and, in consequence, extreme pitching -moment 
corrections may be required. On the other hand, figure 5 which compares the interfer- 
ence factors for the variable y tunnel ^6^ ^ = oj with those for both closed and closed- 
on -bottom -only tunnels (y = 2) indicates that only a relatively small interference would be 
encountered by arbitrarily limiting the maximum width-height ratio to 2.0. An alterna- 
tive procedure is to mount large models slightly below the center of the tunnel. Fig- 
ures 6 and 7 present similar results for the same series of wings when mounted below 
the tunnel centerline by distances equal to 0.1 and 0.2 times the semiheight of the high- 
speed (y = 2) section. Figure 8 shows, as a function of span-width ratio, the required 
distance below the centerline in order to obtain 6 W ^ = 0 at y = 90° and y = 2. It 


11 



■■■mi i hi i ■ 1 1 ii 


can be observed that only small displacements AH/H are required provided that the 
model span is held to moderate values of one-half or less of the wind-tunnel width. 

Effect of model configuration . - Figures 9 to 11 indicate the effect of configuration 
variables on the required schedule of y for minimum interference by comparing the 
computed schedules for several wing sweep angles and for a rotor. For all configura- 
tions, the model span-width ratio a is 0.5 and the angle of attack is zero. For moder- 
ately heavily loaded conditions (x ~ 75°), it can be seen that the configuration considered 
may call for a tunnel width -height ratio as much as 10- to 20-percent different from the 
values given in figures 5 to 7. These differences are calculable, however, and once the 
wind-tunnel procedures are established, it should only be necessary to provide the control 
computer with a table of values computed for the specific configuration being tested. 

Effect of angle of attack .- A change in angle of attack slightly alters the relative 
position of each portion of the model with respect to the tunnel boundaries and conse- 
quently alters the interference caused by these boundaries (ref. 12). In the present con- 
text of a minimal-interference tunnel, this change in interference requires that the 
width-height ratio be altered slightly to retain a zero value of 6 w The magnitude of 
the required changes for a wing (represented by a lifting line) with 45° of sweep and for 
a rotor are shown in figures 12 to 14. Irrespective of mounting height, the changes in 
width -height ratio required by changes in angle of attack are small. Note that for wings 
of lesser sweep than 45°, the effect of angle of attack is even less; for example, an 
unswept wing has correction coefficients which are independent of angle of attack (ref. 12). 

If simplicity of operation is desired, it should be acceptable to simply ignore the 
effect of angle of attack when testing models without tails. On the other hand, it is fea- 
sible, if desired, to supply the computer with a doubly ordered set of width-height ratios 
(as functions of both x and a ) f° r use * n a double interpolation routine. If on-line data 
reduction is anticipated, the latter procedure is preferable, since angle of attack has 
large effects (refs. 12 and 26) on the corrections to pitching moment due to the tail. 

Interference at tail .- Figure 15(a) displays the interference factors in the variable 
width-height ratio tunnel for a zero span model (at a = 0°) at the center of lift l^j H* = 0 
and at tail lengths Z^/H* 0-5, 1-0, and 1.5. These values are compared with the values 

for the equivalent fixed closed and closed -on -bottom -only tunnels having a width-height 
ratio of 2.0. It is evident that, in general, the use of variable geometry decreases the 
interferences by as much as an order of magnitude. It is evident that pitching moments 
due to the tail will be subject to some correction even in the variable -geometry test sec- 
tion; however, the required corrections should be much smaller than in the fixed-geometry 
tunnels. 
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More properly, corrections to pitching moment are nearly proportional to the 
difference in interference factors at the tail and at the center of lift (ref. 26). A similar 
comparison of this difference is presented in figure 15(b). This form of presentation 
does not alter the curves showing the interference at the tail for the variable -geometry 
tunnel; however, ^ is markedly smaller than 6^ ^ for any given tail length in 

the fixed -geometry tunnels since the interference at the center of the lift is of the same 
sign as the interference at the tail in these tunnels. Nevertheless, the variable -geometry 
tunnel still displays remarkably reduced interference effects on pitching moment. Fur- 
thermore, examination of the trend with increasing tail length indicates that the rate of 
change of interference with distance downstream is substantially reduced; thus induced 
camber and other similar effects (ref. 18) will also be reduced. 

Obviously, finite span, configuration differences, and angle of attack may have large 
effects on the absolute values of the results presented in figure 15. Complete discussion 
of these effects would necessitate an extremely protracted discussion. Several sample 
calculations (not presented herein) indicate that the qualitative conclusions already drawn 
are not altered. Readers interested in particular cases are directed to the computer 
programs in the appendixes. 

Distribution of interference factors .- The distribution of interference over models 
in the variable -geometry tunnel is examined and compared with the distributions in fixed- 
geometry tunnels in figures 16 to 19. In each case, the model spans one-half of the tunnel 
width and is centered in the tunnel with a - 0°. The spanwise distribution of interference 
over an unswept wing is shown in figure 16; the spanwise distribution over a wing having 
45° of sweep in figure 17; and the lateral and longitudinal distributions over a rotor in 
figures 18 and 19. In each case, it is obvious that the model in the variable-geometry 
tunnel experiences essentially no interference nonuniformity in dramatic contrast to the 
large nonuniformities which, in general, are found in the fixed -geometry tunnels. Fig- 
ures 16 to 19 collectively indicate that data from the variable -geometry tunnel should be 
essentially free of induced camber and twist effects (ref. 18). 

Re sidual interf er enc e . - It is noted that the variable -geometry feature of the test 
section has been used to reduce the vertical interference due to lift (6 W to zero. 

There remains a vertical interference due to drag (5 w jj) and streamwise or horizontal 
interferences due to both lift (6 U and drag (6 U . In general, these interferences do 

not become zero simultaneously with 6 W however, the effect of these interference 
factors on the data is so much less than the effect of 6,„ T that their nonzero values 

W, J-j 

are not too troublesome. 

Figure 20 compares the values of the residual interference factors in the variable 
and fixed-geometry tunnels. Except for conditions near y = 90°, where these factors 
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are essentially zero in the fixed -geometry closed tunnel, the values are somewhat smaller 
in the variable -geometry tunnel. The most significant reduction that results from the use 
of variable geometry is in 6* ^ at low speeds (low y); under such circumstances, this 
interference factor can result in significant reductions in the effective speed of a fixed 
tunnel. 

If the induced drag-lift ratio of the model is known within reasonable limits prior 
to the test, it is possible to develop a schedule of width -height-ratio that reduces all 
vertical interference to zero at the model, leaving only a velocity (or "q") correction. 

The total vertical interference velocity is (from refs. 3 and 4) 


Aw = 


J w,L 


“m 

A'p 


w 0 + 6 w,D 


•"■m 

Aq 1 


U r 


However, from reference 27 



so that equation (2) may be written as 


Aw = 


D i 

6 w,l + r 5 


\ A m 

w,dJa7 W ( 


( 2 ) 


(3) 


(4) 


For a rotor at low or moderate tip-speed ratios, where the resultant force is essentially 
perpendicular to the rotor disk, equation (4) may be written as 

Aw = ( 8 w,L + 6 w,D tan w 0 < 5 > 

Provided only that (or, in the case of a rotor, a) is known within reason- 

able limits, it is possible to modify the programs in the appendixes to reduce the 
bracketed sum in parentheses in equation (4) or (5) rather than just 6 w L to zero. This 
refinement has not been provided in the programs provided herein; however, the required 
modifications are straightforward. 

Effect on recirculation li mits.- The most significant recent development in low- 
speed V/STOL wind-tunnel testing has been the discovery (refs. 23 and 24) of minimum 
speed limits below which the flow in the tunnel recirculates in a manner unrepresentative 
of free air. These limits can be expressed (ref. 18) in terms of a minimum distance 
behind the model for the intersection of the floor and the theoretically straight wake. The 


14 



exact minimum distance depends upon the wind-tunnel proportions and certainly upon the 
tunnel configuration (for example, the presence of corner fillets has a strong influence). 

The investigations of references 23 and 24 explored recirculation limits for closed 
tunnels only. It appears, however, that some similar restriction will hold in almost any 
wind tunnel even though the value of the limiting impingement distance is not known at 
present except for the closed tunnels. In this paper, discussion of recirculation limits 
is confined simply to the ratio of the impingement distance in the variable -geometry sec- 
tion to that in a fixed -geometry section. Even though such a simplistic presentation 
ignores many known effects, it at least provides a qualitative guide to the merit of a par- 
ticular arrangement with regard to minimum test speeds. 

For a model mounted in the center of the tunnel, the impingement distance behind 
the model is 


Xf = H tan x (6) 

In the present case, where the width of the test section is fixed, H = B/y; so that 

_ ^ tan X _ r* _ 2 (7) 

*f* ^tan x y y 

Figure 21 shows the variation of x^/x^* as a function of wake skew angle for sev- 
eral unswept wings of differing span-width ratio. The schedules of y previously given 
by figure 5(a) are used. The large-span wings of this group are shown to be at a dis- 
advantage with respect to recirculation at the higher wake skew angles. As noted pre- 
viously, however, only a very mild correction is encountered if the width-height ratio is 
limited to a maximum of 2.0. If such a limit is imposed, the curves of figure 21 would 
be at Xf/xf* = 1.0 until the curves shown actually cross that line. When minimum speed 
is the primary concern, this procedure is preferable to the alternate procedure of 
mounting the model below the centerline since a low mounting height also reduces xj. 
Except for the high -skew -angle region just discussed, it is obvious that the use of a var- 
iable width -height-ratio tunnel greatly increases the impingement distance ratio - by 
factors as great as 4 — and thus allows recirculation-free testing to substantially lower 
speeds than are allowable in a fixed-geometry tunnel. 

Additional applications of variable y tunnel .- The same features required to pro- 
vide the variable geometry in the tunnel can be used to add several additional test capa- 
bilities. A few of these capabilities are sketched in figure 22. 
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Asymmetric motion of the entrance cone and diffuser can provide an inclined-flow 
facility suitable for unpowered free-flight tests of models with low lift-drag ratios, such 
as parawings and lifting reentry bodies. No equivalent large-scale test facility exists at 
present and such work must normally be conducted under less controlled conditions in 
actual flight tests. 

Both entrance cone and diffuser can be moved upward simultaneously to provide for 
ground-effect testing. This procedure is far preferable to lowering the model if the wind 
tunnel employs an external mechanical balance system since it obviates the need for 
moment transfers through the large distances between the balance and model moment- 
centers. 

Either the floor alone or the floor together with the entrance and exits can be raised 
to the level of the model to provide a convenient height at which to work on the model. 

The value of this feature in a very large tunnel will be apparent to anyone who has ever 
worked in one of the so-called full-scale tunnels. 

Closed -On -Bottom -Only Tunnel With Variable Model Height 

General concept of configuratio n.- The second type of variable -geometry wind- 
tunnel test section to be discussed is the closed-on -bottom -only tunnel with fixed bound- 
aries but with the model mounting height programed as a function of wake deflection. 

This configuration is of particular interest since a tunnel utilizing a test section of this 
type having a width-height ratio of 1.5 is presently under construction at the NASA Langley 
Research Center (ref. 28). Figure 23 is a photograph of a model of this test section, 
illustrating not only the general type of test section under consideration herein, but also 
one manner in which removable walls may be employed in order to provide complementary 
closed or slotted sections for othe^ypes of work. 

The variable model-height tunnel is considerably simpler and cheaper to construct 
than the variable width-height ratio tunnel previously discussed. In many tunnels 
employing sting balances, most of the mechanical elements required for this class of 
tunnel are already present since sting mounts are generally capable of considerable ver- 
tical motion. Essentially the only requirements (other than removable walls and ceiling 
and an exit collector cone) are the provision of a control computer with suitable velocity 
and force transducers. No intrinsic design problems appear to exist other than insuring 
a stable air stream without severe oscillations or pulsations. In this regard, there is 
considerable past work (for example, refs. 29 and 30) to provide, guidance in obtaining a 
clean basic test-section flow. 

Computer programs .- Computer programs for the variable -model -height tunnel are 
also included in the appendixes. Appendix C provides a computer program for wings; 
appendix D provides a similar program for rotors. 
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Required schedule of model height for wings .- Since the wall effects and their 
changes with vertical model position are both functions of the width-height ratio of the 
wind tunnel, it is necessary to examine, at least initially, a family of tunnels with variable 
model height. Figures 24, 25, 26, and 27 show the required schedule of model height for 
unswept wings of varying span -width ratio in tunnels having width-height ratios of 2.0, 

1.5, 1.0, and 0.5, respectively. In each case, the resulting interference (6 w ^ = 0) is 
compared with the interference factors in both closed and closed-on-bottom-only tunnels 
of the same width-height ratio when the model is retained in a central mounting position. 

For the cases studied, the maximum range of vertical motion is not extraordinarily 
large, ranging from about 0.6 of the tunnel semiheight H for the widest (y = 2.0) tunnel, 
to about 0.8 H for the deepest (y = 1.0 and 0.5) tunnels. The required positions are, how- 
ever, lower in the tunnel as the width-height ratio decreases. In the most extreme case 
(y = 0.5 and cr = 0.75), the model would be within one-quarter of a semiheight of the floor 
in the high-speed (y = 90°) condition. This result should be expected since the deeper 
tunnels of this type have a greater ratio of open-boundary perimeter to closed-boundary 
perimeter than the wider tunnels; consequently, the model must be moved closer to the 
solid boundary in order to increase effect of the floor (upwash) to the point where it is 
exactly equal to the effect of the open boundaries (downwash). 

The consequences of the proximity to the floor in the deeper tunnels will be explored 
more thoroughly in succeeding sections of this paper. It is obvious at the outset, how- 
ever, that the deep tunnels of this class will present certain operational problems with 
regard to clearance of (and interference at) the tail with large angles of attack as well as 
problems in reaching low speeds without violating recirculation limits. 

Interference at tail .- The effect of variable model-height operation on the interfer- 
ence at the tail in the same family of tunnels is explored in figures 28 to 31. Compared 
to a fixed central location, the effectiveness of variable model height in reducing the cor- 
rections to pitching moment decreases with both decreasing width -height ratio and 
increasing tail length. Indeed, in the deepest tunnel (y = 0.5), the correction to pitching 
moment caused by the tail would be decidedly increased, and in the square tunnel (y = 1.0), 
on balance, the same order of correction would be obtained. 

Variable model -height operation does not reduce pitching-moment corrections to the 
same low level as could be achieved in the closed-on-bottom-only tunnel of variable width- 
height ratio (fig. 15). This relative disadvantage is a price that must be paid for the 
greater simplicity and lesser costs of the variable model-height configuration. 

In the wider tunnels of this family (y = 2.0 and 1.5), it is observed that a significant 
improvement in pitching -moment interference is obtained provided that the tail length is 
not extreme. Modest tail lengths on the order of l\/H = 0.5 will require pitching- 
moment corrections of about one-half those normally encountered. 
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The fact that the improvement in A5 W L is greatest for the shortest tail lengths 
is significant since it indicates that the longitudinal gradient of interference velocity at 
the model lift-center is probably greatly reduced. Thus, this type of tunnel should pro- 
duce only small induced-camber effects. 

Residual interference .- The effect of variable model-height operation on the resid- 
ual interference factors <5„ T , 6 and 6 r, for the same four tunnels is shown in 

U,1j U. y W jU 

figures 32 to 35. No truly significant change in the level of these factors is evident. As 
in the case of the variable width-height ratio tunnel, the programs given in appendixes C 
and D can be modified to eliminate the vertical interference due to drag provided that the 
drag-lift ratio of the model can be approximated within reasonable limits. 

Effect of model configuration .- The effect of model configuration as exemplified by 
a rotor and wings of varying sweep, all having span -width ratios of 0.5, are illustrated in 
figure 36 for one width-height ratio (y = 1.5). (Results for the remaining three width- 
height ratios are omitted for conciseness in this and the following several sections. 
Equivalent information for other width -height ratios can be obtained, if required, by use 
of the programs given in appendixes C and D.) At least for y = 1.5, figure 36 indicates 
that the effect of configuration is so small that it could be neglected unless extreme pre- 
cision is required in the data. 

Effect of angle of attack . - The effect of angle of attack on the required model- 
height schedule for minimal interference is illustrated in figure 37. The same tunnel 
and span are considered as in figure 36. Figure 37(a) presents the results for a wing 
with 45° of sweep; figure 37(b) presents the results for a rotor. In the cases considered 
herein, it is obvious that the effect of angle of attack is negligible. On the other hand, it 
has been observed previously that this type of variable -geometry tunnel will require cor- 
rections to pitching moment. As noted in references 12 and 26, angle of attack would be 
expected to have significant effects on this pitching-moment correction. 

Distribution of interference .- The effect of variable model-height operation of the 
closed-on-bottom-only tunnel on the distribution of interference for wings having 0° 
and 45° of sweep is illustrated in figures 38 and 39. Corresponding results for the dis- 
tribution of interference along the lateral and longitudinal axes of a rotor are given in fig- 
ures 40 and 41. It will be noted in each case that the interference distribution is usually 
better in the variable model -height tunnels than in the tunnels with fixed, centrally located, 
models. In particular, the longitudinal distribution of interference over the rotor (fig. 41) 
is significantly flattened by properly varying the model height. This result is in accor- 
dance with the previous observations about the longitudinal gradient of interference dis- 
cussed earlier in relation to the interference at the tail. 
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Effect of recirculation limits .- The comparison of the variable model -height closed- 
on-bottom -only tunnels of different width-height ratios with respect to relative recircu- 
lation limits is more complex than might be apparent at first glance. Note that in this 
case, the impingement distance is 


Xj = h tan x 


( 8 ) 


and, since £ = H/h 

H tan x 

*f _ g 1 _ AH [ x 

x f) t H tan x ? H 
1 


Equation (9) may be used to compare the effect of using variable model height in 
tunnels of any given width-height ratio. The results of such a comparison are presented 
in figure 42 for unswept wings of varying span in tunnels having each of the four different 
width-height ratios considered herein. In general, it would appear that the use of a 
variable model height is deleterious at the high skew angles but that it may be heipful in 
the wider tunnels if the wake is depressed sufficiently from the horizontal. 


In comparing the various width -height ratios, it should be noted that if all of the 
tunnels have the same area (which, for a given maximum speed, generally determines 
the cost of constructing the tunnel), the tunnel semiheight (used in eq. (9)) will be different 
in each tunnel. Note that the tunnel cross-sectional area is = 4BH = 4yH^ so that 
the ratio of semiheights for equal is 



( 10 ) 


Thus, proceeding as before, the following result is found for equal area tunnels 






an 


A comparison on this basis is presented in figure 43 for models of fixed span-width 
ratios of 0 and 0.5. Overall, it would appear from figure 43 that the square tunnel 
(y = 1.0) of this class is superior from a viewpoint of recirculation limits. More 
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realistically, however, an experimenter designing a model for testing will find that his 
model size is physically fixed and that it is not possible to "rubberize” the model design 
to obtain a fixed span-width ratio. In general, the physical size of the model is fixed 
rather than the span -width ratio. 

Observe that the model semispan in tunnels of constant is 

s = oB = oyH = oy(H/H*)H* = ftyH. *a 


and 


s \f2ya 
s * ~ \[ 2 { 2)o* 


or 



Thus, for fixed span models in fixed area tunnels, both H/H* and cr/a* must 
vary in the identical manner (eqs. (10) and (12)). On this basis, the comparison of the 
various tunnels is unaltered from figure 43(a) when a = cr* = 0; however, for finite span- 
width ratios of cr* = 0.25 and 0.5, the comparison of the different tunnels is as in fig- 
ure 44. Note that the result for y = 0.5 and a* = 0.5 is not shown since in that case 
the model wing tips actually touch the side boundaries. Examination of figure 44 indicates 
that, on an overall basis, the wide tunnels of this class are probably superior from the 
viewpoint of recirculation limits for models of practical spans. The penalties of variable 
model -height operation in these tunnels are very small and are confined to the higher 
speed conditions (high y) where recirculation is not likely to be a limiting factor. 

For severe conditions (low y), the wider tunnels provide a substantial improvement 
in the relative recirculation limits. On the other hand, it is observed that the improve- 
ment found in this case is significantly less than that obtained in the variable width-height 
ratio tunnel discussed previously (fig. 21). 

Choice of width-height ratio . - In view of the conflicting requirements of different 
types of wind-tunnel testing, particularly when nonaeronautical applications (that is, 
forces on bridges, trains, etc.) are considered, the choice of the optimum proportions for 
a wind tunnel is always somewhat imprecise and intuitive. In the present paper, where 
very low-speed tests of V/STOL vehicles are the primary concern, the impact of recir- 
culation limits on tunnel utility and useful speed-range would appear to rule out the choice 
of a deep tunnel. The knowledge that pitching -moment corrections due to interference 
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at the tail are not reduced in the square tunnel by variable model -height operation makes 
the square tunnel significantly less attractive. Furthermore, the fact that a model of 
given span has a smaller span -width ratio when tested in a tunnel of larger width -height 
ratio (thus reducing the lateral nonuniformity of wall interference) indicates that the 
tunnel should be substantially wider than square. 

If the width-height ratio is chosen to be very large, then, for fixed tunnel areas, the 
floor is found to be quite close to the model. In a practical sense, the clearance from 
the floor of the model tail at high angles of attack can be restrictive. In addition to the 
problem of physical clearance, the result is interference at. the tail which is a severe 
function of angle of attack; in certain cases such interference has been found to produce 
increases in apparent longitudinal static stability that were of the same order as the 
basic stability of the model. 

On an overall basis, it would appear that a width -height ratio on the order of 1.5 
would be satisfactory and would still yield most of the advantages of variable model -height 
operation. This is the value chosen for the tunnel presently being constructed at the 
Langley Research Center. It is noted in passing that width-height ratios substantially 
equal to 1.5 (that is, 7 x 10 and 8 x 12) have generally proven quite satisfactory in the 
past. 

Operation for zero moment-correction . - Since one of the main shortcomings of the 
variable model -height configuration is the minor relief from pitching -moment correction, 
it is interesting to examine the possibility of operating the tunnel so as to minimize 
pitching-moment corrections rather than the corrections at the center of lift. Figure 45 
shows the interference factors 6 W ^ for both the wing and the tail for one tail length 
and one angle of attack as calculated by the method of reference 12. The circled points 
in figure 45 show the intersections at which the two interferences are equal so that 
A6 W ^ = 0. For these conditions, there is essentially no correction required for pitching 
moment due to the tail. 

This procedure has been carried through graphically (but in more detail) for a 
variety of width-height ratios, tail lengths, and angles of attack. The results are pre- 
sented in figures 46 to 49. The most significant observations to be drawn from these fig- 
ures is that operation for zero moment-correction is only possible for a limited range of 
skew angles, tail lengths, and angles of attack; the range depends to a large extent upon 
the width-height ratio, being largest for the tunnels of greatest depth. Furthermore, the 
corrections at the center of lift can become very large in the downwash (6 W ^ > Oj direc- 
tion. The corrections at the center of lift, in general, are significantly larger in absolute 
magnitude than when operating with the model fixed at the center of the tunnel. (Compare 
figs. 46 to 49 with figs. 24 to 27.) 
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Except under unusual circumstances, it would appear that operation in this mode is 
contraindicated. Consequently, no further examination of this possibility is included 
herein, and no computer programs similar to those in the appendixes have been developed. 

Closed Tunnel With Variable Width -Height Ratio 

General concept of configuration .- Even though it is not possible to reduce S w ^ 
to zero, it is of value to explore the extent to which interference can be reduced by 
applying variable-geometry concepts to closed tunnels. The first type of closed tunnel 
to be discussed is the variable width -height-ratio closed tunnel. 

In appearance, such a tunnel would look very much like that shown in figure 3 
except that the walls of the two-dimensional contraction would continue along the sides of 
the test section and a closed ceiling would be carried with the upper lips of the entrance 
and exit. The operational features, additional applications, and developmental problems 
of such a tunnel are essentially the same as those discussed previously for the similar 
closed-on-bottom-only tunnel. One advantage to the closed version is that it is consid- 
erably simpler to achieve a clean nonpulsating flow in a completely closed test section. 
Furthermore, the boundary conditions to be satisfied at all the boundaries are more 
firmly known, since there are no small perturbation assumptions in the development of 
these conditions (ref. 31). As a partial balance against the firmness of the boundary con- 
ditions, the possibility of separation from the ceiling under conditions of extreme force 
coefficients (ref. 1) should be pointed out. 

Computer programs for variable y closed tunnel .- Only minor modifications 
(largely to reduce input values to correspond always to the high-speed-section values) 
are required to use the programs developed in reference 12 for calculation of the inter- 
ferences in the tunnel under discussion. Appendixes E and F present programs for wings 
and rotors with these modifications already made. 

Optimum schedule of y.- The interference factor 6* T is presented in figure 50 
for several different models and for a range of width-height ratios. This presentation is 
strikingly different than that of figure 1(a) since the increase in tunnel .area with 
decreasing width-height ratio is accounted for by the use of 6* T rather than 6 T . 

W y J_l W y i-J 

Since the interference decreases monotonically with y, it is obvious that the opti- 
mum manner of operation is to decrease y (thus increasing cross-sectional area) as 
rapidly as the required test speed permits. Thus, all speed control in the upper speed 
range of the tunnel could be vested in the control of the tunnel area rather than in either 
fan rotational speed or fan pitch. This feature could significantly simplify design of the 
speed control system in the tunnel. 
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As was indicated in figure 1(a), the interference factors 6 w ^ increase after the 
width-height ratio is decreased below some optimum value. It is obvious from figure 50 
that the increase in 6 w L becomes big enough to largely overcome the advantage of the 
increases in area for width-height ratios on the order of 1.0 or less; that is, the reduc- 
tion in 6^ ^ is small for reducing y to values much less than 1.0. Thus, it would 
appear from figure 50 that the tunnel need only have a range of y from 1.0 to 2.0 in 
order to achieve essentially all of the possible beneficial results of this type of operation. 
Other factors modify this conclusion to a certain extent as will appear in subsequent 
discussion. 

Interference for particular models . - It is obvious from examination of figure 50 
that the interference experienced by a given model will be dependent upon its own varia- 
tion of skew angle with forward velocity as the tunnel area and speed change with changes 
in width-height ratio. If it is assumed for simplicity that the power required is propor- 
tional to ApV 2 (equivalent to assuming that the energy ratio is unaltered by changes in 
the tunnel opening), and A-p = 4BH = (4/y)B 2 , so that for a constant-power opening of the 
tunnel 


Power ~ — B 2 (V*) 3 = - B 2 V 3 

v* y 


(13) 


Solving equation (13) for V, nondimensionalizing the result by w^, and taking 
y* = 2 yields 


_v_= Cl ) 1 / 3 Xl 

w h \ 2 / w h 


(14) 


From reference 27, with D^/L taken as zero (note that the assumption that = 0 
may require the presence of some negative Dj, or thrust, from some portion of the 
model) 



and, also from reference 27 



cos y, 


m 


(16) 
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Substitute (16) into (15) and solve for cos x m to obtain 



Then the "effective skew angle" (ref. 9) may be taken with sufficient accuracy as 


*m + 90 ° 


The foregoing result is satisfactory for "aerodynamic systems" in which the forces 
are produced by changes in velocity and direction of a constant mass of air which enters 
at infinity upstream and leaves at infinity downstream. These conditions are violated in 
the case of cold jets supplied from within a model to simulate jet engine flows. For such 
systems, equation (14) is still valid; however, the momentum skew angle is more appro- 
priately taken from the vector diagram in the sketch. 




( 20 ) 
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Actually, of course, the wake of the jet must leave the model tangent to the exit 
nozzle and then curve rearward until the inclination to the free stream is negligible at 
very large distances downstream. The results of reference 18, as well as a comparison 
with the flow photographs of reference 32, indicate that the appropriate effective skew 
angle is again as given by equation (19). 

The relationship between x ^ 7 is given by the solution of equation (18) (or 
eq. (20), if appropriate) and of equation (19) and depends only upon the reference veloc- 
ity wjp which in turn (ref. 27) is only a function of the equivalent "disk-loading." The 
dashed lines in figure 50 represent the paths followed by models having several differ- 
ent ratios of w^/V*. For a 300-knot tunnel, the curve labeled w^/v* = -0.06 might 
represent either a helicopter rotor or a very heavily flapped wing; the curve labeled 
wjj/v* = -0.2 might represent a tilt-wing VTOL aircraft; the curve labeled w^/v* = -0.6 
might represent a fan-in-wing design; and the curve labeled wj^V* = -2.0 might repre- 
sent a jet-lift VTOL design. (The first three cases were computed using equation (18); 
the final case used equation (20).) 

In the highest speed condition, the tunnel is always operating under the same inter- 
ference conditions that it would experience if the tunnel were a fixed -geometry tunnel 
with y = 2; however, as the speed is reduced, the interference reduces rapidly to values 
that are on the order of one-half to one-third those of the fixed tunnel. The entire opening 
from y = 2.0 to y = 0.5 occurs in the upper one-third of the speed range of the tunnel. 
For all lower speeds, the interference is given by the curves for y = 0.5. It is noted 
that during operation at points above the maximum opening, the tunnel velocity is com- 
paratively large. Under these conditions the force coefficients (and, consequently, uq 
and wq) are small so that the interference of the walls is small. It is at low speeds, 
where the force coefficients and mean induced velocities are large, that wall effects 
become truly significant. It is precisely this range in which the variable width-height- 
ratio operation of the closed tunnel reduces the wall effects. 

Interference at the tail . - The interference at the tail, in terms of A 5^ L , is shown 
for zero span and for three tail lengths in figure 51. It will be noted that, for low speeds, 
corrections for pitching moment at the tail will be extremely small as a result of 
variable -geometry operation. In particular, if the tunnel opens to as small a width-height 
ratio as 0.5, such corrections almost vanish; this tunnel reduces tail corrections to even 
lower values than the variable y closed-on-bottom-only tunnel previously discussed. 
(Compare fig. 15 and fig. 51.) The greater reduction in the present case is entirely due 
to the more rapid opening schedule. 

Residual interference .- Figure 52 presents the effect of variable-geometry opera- 
tion of this tunnel on the interference factors 6* L , 6* D , and 6^ q. It is evident that 
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these factors are also drastically reduced. In particular, it is noted that the streamwise 
interference is reduced to negligible values at low speeds. 

Distribution of interference.- Figures 53 to 56 present the distribution of 6* T 

— — 1 iV j Xj 

over the same models studied with respect to the previously discussed wind-tunnel con- 
figurations. If the curves for y - 2.0 are compared with the curves for the appro- 
priate y (taken from the intersections of the solid and dashed curves of fig. 50), it is 
obvious that variable y operation of the closed tunnel results in dramatic reduction of 
the nonuniformity of interference at low speeds. 

Effect on recirculation limits .- The effect of variable y operation on the impinge- 
ment distance is shown in figure 57 as a function of either y or V/V*. It is possible, 
of course, to include in this figure the effect of y on the allowable recirculation distance 
(refs. 18 and 23); however, this has not been included in figure 57 primarily in order to 
provide a more consistent comparison with the previously discussed tunnel configurations. 

The rapid opening schedule of the variable y closed tunnel results in sudden and 
rapid increases in impingement-distance ratio. The extent of the increase is limited 
only by the maximum opening of the tunnel. If the tunnel is intended primarily for work 
in the low -speed V/STOL transition regime of flight, it is obvious that the maximum 
desirable opening (or lowest y) is the maximum opening that is economically feasible. 
When the minimum value of y is about 2/3 or less, it is obvious that this type of tunnel 
will provide tests free of recirculation effects at speeds (or y) lower than any other 
tunnel considered herein. 

Closed Tunnel With Variable Model Height 

It remains to consider the possibility of variable model -height operation of a closed 
tunnel. For each of four width -height ratios, the schedule of model height required to 
produce minimum 6,„ T (determined for a = 0 from the computer programs developed 
in ref. 12), is shown, together with the effect on 6 w ^ in figures 58 to 61. It is evident 
that the reductions in interference achieved by this technique in the closed tunnel are only 
minor and are probably not worth the expense and difficulty of achieving the required 
motion of the model. 


CONCLUSIONS 

This study of variable -geometry wind tunnels, using changes in either width-height 
ratio or model height, to minimize boundary interference in wind-tunnel tests of V/STOL 
aircraft has indicated that several promising possibilities exist: 

1. The closed-on -bottom -only configuration with properly scheduled width-height 
ratio can reduce the vertical interference velocities to zero at the lifting system itself. 
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Simultaneous major improvements in uniformity of interference, pitching-moment correc- 
tions, and minimum speed for recirculation-free testing are obtained. 

2. Closed-on-bottom-only tunnels, wider than square, can also be used to reduce 
the vertical interference velocity at the lifting system to zero if the model height is 
properly scheduled. Improvements in uniformity, pitching -moment corrections, and 
minimum recirculation-free test speed are also obtained; however, these latter improve- 
ments are significantly less than those obtained in the variable width -height ratio 
configuration. 

3. At low speeds, closed tunnels having variable width-height ratio offer the possi- 
bility of reducing wall effects by factors of 2 or 3 at the lifting system. Such tunnels 
reduce, at low speeds, the nonuniformity of interference and pitching -moment correc- 
tions to the vanishing point. These tunnels can also provide the lowest speeds for 
recirculation- free testing. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., April 24, 1969, 

721-01-00-20-23. 
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APPENDIX A 


FORTRAN PROGRAM FOR COMPUTING INTERFERENCE OVER WINGS 
IN A VARIABLE WI DTH-HE I GHT-RAT 10 CL OS EC-ON- BO T TOM -ON L Y TUNNEL 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN, VERSION 2.1, TO RUN ON CDC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


NOTE THAT ALL INPUT AND OUTPUT VARIABLES ARE REFERENCED TO THE GAMMA=2*0 
TEST SECTION. INPUT IS AT ADDRESS 1 IN FORMAT 900 (TWO CAROS PER CASE). THE 
REQUIRED INPUT VARIABLES ARE 


LI SPAN LOAD-DISTRIBUTION INDICATOR, LI=1 FOR UNIFORM LOADING, 

L I -2 FOR ELLIPTIC LOADING 

ZETA1 SEMIHEIGHT OF GAMMA =2 SECTION OIVIDED BY HEIGHT OF AERODYNAM I C 

CENTER ABOVE FLOOR IN GAMMA=2 SECTION 

ETA1 DISTANCE FROM AERODYNAMIC CENTER TO RIGHT-HAND WALL DIVIOED 

BY TUNNEL SEMIWIDTH 

SIGMA RATIO OF WING SPAN TO TUNNEL WIDTH 

LAMBDA WING SWEEP ANGLE, DEG 

ALPHA ANGLE OF ATTACK OF WING, DEG 

ITAIL NUM8ER (0 TO 9) OF TAIL POSITIONS, PROPORTIONATELY INCREASED 

IN TAIL LENGTH AND HEIGHT, AT WHICH INTERFERENCE IS 
REQUIRED. USE SMALLEST LENGTH AND HEIGHT AS INPUT. 

SIGMAT RATIO OF TAIL SPAN TO TUNNEL WIDTH 

TL DISTANCE OF TAIL BEHIND AERODYNAMIC CENTER, NONDI ME NS IONAL I ZED 

WITH RESPECT TO SEMIHEIGHT OF GAMMA=2 SECTION 

TH HEIGHT OF TAIL ABOVE AERCDYNAMI C CENTER, NONDI MENS I ONAL IZ ED 

WITH RESPECT TO SEMIHEIGHT OF GAMM A=2 SECTION 


THIS PROGRAM REQUIRES SUBROUTINE DLTAS (SEE APPENDIX 
VARLET (SEE APPENDIX H). 


G) AND SUBROUTINE 


PROGRAM W I NDTUN ( INPUT, OUTPUT ,TAPE5=I NPUT ,TAPE6=0UTPUT ) 
COMMON ZETA, ETA, GAM MA , XOVERH , YOVERH , ZOVE RH , DE LTA ( 28 ) 
DIMENSION XLEilO ) »XLQAD( 10) ,XDELTM28 ) »C( 8J,ST0RE(8) 
REAL LAMBDA 


(A 1) 
(A 2) 
(A 3) 
(A 4) 
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DATA (C( I) » I s 1 v 8 )/ 20., 30 .,40., 50., 60., 70., 80., 90./ 

XL E ( 1)=XLE( 10>=0. 43579 
XLE(2 )=XLE(9) =0.71422 
XLE( 3)— XLE (8 ) =0 • 86603 
XLE(4)=XLE(7 )=C • 953 94 
XLE(5)=XLE(6)=0. 99499 

l READ (5,900) LI, ZETA1, ET A1 ,S IGMA , LAMBDA, ALPHA,ITAIL, SIGMAT ,TL,TH 
IF (EOF, 5) 999,47 
47 IF (LI.EQ.l) GO TO 804 
IALPHA=8HELLIPTIC 
SUML=0. 0126104 
DO 808 M2 = 1,10 

808 XLOAD(M2) = XLE(M2 ) 

GO TO 160 

804 S(JML = 0.01 

I ALP HA=8HUN I FORM 
DO 809 M2= 1 , 10 

809 XLOAD ( M2 ) = 1 * 0 

160 WRITE (6,901) SIGMA, IALPHA, LAMBDA, ZETA1, ETA1, ALPHA 
WRITE (6,210) 

WRITE (6,211) 

WRITE (6,212) 

WRITE (6,213) 

WRITE (6,214) 

WRITE (6,215) 

WRITE (6,216) 

WRITE (6,217) 

WRITE (6,218) 

AALP= ALPHA 

AL AM= LAMBDA 

AT L = TL 

ATH=TH 

AZ ET A = ZE T A 1 

RAD=0.0174532925199 

ALPHA= ALPH A*RAD 

L AMBDA=L AM BDA*R A D 

CONST 1 = 1 .0 

DO 41 K = 1,8 

DO 803 L 1= 1 , 28 

DELTA(L1 ) =0 • 

803 XDELT A( LI ) =0. 

GAMMA=GMA1=2.0 
I NK= 1 

IF ( SIGMA. NE.O. ) GO TO 811 
M7=N7 = 1 
XLOAD ( 1 ) = 1 .0 
SUML = 1 • 

GO TO 812 

811 IF (ETA1.NE.1.) GO TO 813 
M7 = 5 

N7 = 10 
CONST 1 = 2 • 

GO TO 812 

813 M 7=N 7= 1 0 

812 AC=0.5*S IGMA*TAN( LAMBDA) *GAMMA 
IF (LI.NE.l) AC =0 • 84 8*A C 

ZETA1 = AZETA/ (1.0 +AC *AZETA*SIN(ALPHA ) ) 

ZETAl=2.0/(2*0+GAMMA*( 1.0-ZETA1) ) 

DO 801 Ml= 1 , M7 

DO 802 N 1 = 1 , N7 

XSTAR=( 11.-2. *FL OAT (Ml) ) /10. 


(A 5) 
(A 6) 
(A 7) 
(A 8) 
(A 9) 
(A 10) 
(A 11) 
(A 12) 
(A 13) 
(A 14) 
(A 15) 
(A 16) 
(A 17) 
(A 18) 
(A 19) 
(A 20) 
(A 21) 
(A 22) 
(A 23) 
(A 24) 
(A 25) 
(A 26) 
(A 27) 
(A 28) 
(A 29) 
(A 30) 
( A 31 ) 
(A 32) 
(A 33) 
(A 34) 
(A 35) 
(A 36) 
(A 37) 
(A 38) 
(A 39) 
(A 40) 
(A 41) 
(A 42) 
(A 43) 
(A 44) 
(A 45) 
(A 46) 
(A 47) 
(A 48) 
(A 49) 
(A 50) 
(A 51) 
(A 52) 
(A 53) 
(A 54) 
(A 55) 
(A 56) 
(A 57) 
(A 58) 
(A 59) 
(A 60) 
(A 61) 
(A 62) 
(A 63) 
(A 64) 
(A 65) 
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YSTAR=(2.*FL0AT( N1 1-11.)/ 10. 
ZSTAR=(ll.-2.*FLGAT(Ni))/10. 

ET A=ETA1+YSTAR*S IGMA 

ZETA=ZETAl/( l.-A8S( YSTAR ) *$ I GMA*G AMMA+ZE TA1 *T 
1 ) 

XOVERH=SlGMA*GAMMA*TAN(LAMBDA)*CCS( ALPHA )*( AB 
Y0VERH=f FLOAT (Ml ) -FLOAT! N1 ) ) *S IGMA*GAMMA * ( -. 2 
ZOVERH=SIGMA*GAMMA*TAN(LAMBDA)*S IN ( AL PHA ) * ( AB 
IF ( INK. EQ .7 ) GO TO 850 
CALL VARLET (C< K ) ) 

********************* see APPENDIX H FOR SUBROUTINE 

GO TO 832 

850 CALL DLTAS ICC K I I 

********************* SEE APPENDIX G FCR SUBROUTINE 

832 DO 805 Ll=l,28 

805 XDELTA(Ll)=XDELTA(Ll)+( DELTA! LI )*XL0AD(N1) ) 
802 CONTINUE 
801 CONTINUE 

IF < INK.EQ.7 ) GO TO 835 
DELTA(5)=XDELTA( 5 I *SUML*CONS T 1 
IF (ABS! DELTA15) ) .LT. 0.00005) GO TO 833 
DEL2=DELTA ( 5 ) 

IF ( INK.NE.l) GO TO 831 
DELl=OELTA( 5 ) 

GAMMA=GMA2=1 .0 
I NK=2 

XDELT A I 5 ) = 0 • 0 
GO TO 812 

831 IF ( ABS(DELTA(5) ).LT. 0.00005) GC TO 833 
GAMMA=GAMMA- ( GMA 2-GMA1 ) *DEL2 / ( DEL2-DEL 1 1 
IF (GAMMA. LE. 0. 1 ) GAMMA=0. 1 
DEL1=DEL2 
GMA1=GMA2 
I NK=3 

GMA2=GAMMA 

IF (GMA2.EQ.GMA1 ) GO TO 836 
XDELT A ( 5 ) = 0 • 0 
GO TO 812 
836 STORE ( K. ) -0 .01 
GO TO 41 

833 I NK= 7 

XDELT A( 5 ) =0.0 
STORE(K)=GAMMA 
GO TO 812 

835 00 807 L 3= l » 28 

807 DE LTA ( L3 ) = XDELT A (L3 )*SUML*C0NSTl*(GAMMA/2. ) 
WRITE (6,149) C ( K ) , GAMMA 
WRITE (6,150) (DELTA(I), 1 = 1,25,4) 

WRITE (6,151) (DELTA(I), 1=2,26,4) 

WRITE (6,152) (DELTA( I ), 1=3,27,4) 

WRITE (6,153) (DELTA! I ), 1 = 4,28,4) 

IF ( SIGMA. EQ.O. ) GO TO 822 
IF (C(K).EQ.90.) GO TO 820 

XSl=TAN(C(K)*RAD)/( 2 .0*S I GMA*GAMMA*Z ETA1 ) 
CHIM=2.0*C(K)— 90. 

XS2=TAN(CHIM*RAD )/ ( 2 . 0*S I GMA*GAMM A*Z ETA1 ) 

WRITE (6,120) XS1,XS2 
GO TO 822 
820 WRITE (6,119) 

822 DO 814 L4=l , 28 
814 XDELTA ( L4 ) =0 . 


(A 66) 
(A 67) 
(A 68) 

AN(LAMBDA)*SIN( ALPHA) (A 69) 

(A 70) 

S ( X STAR )— ABS ( ZSTAR ) ) (A 71) 

) (A 72) 

S ( Z STAR )—AB S ( X$ TAR ) ) (A 73) 

(A 74) 
(A 75) 

VARLET ********************* 

(A 76) 
(A 77) 

DLTAS ********************** 

(A 78) 
(A 79) 
(A 80) 
(A 81) 
(A 82) 
(A 83) 
(A 84) 
(A 85) 
(A 86) 
(A 87) 
(A 88) 
(A 89) 
(A 90) 
(A 91) 
(A 92) 
(A 93) 
(A 94) 
(A 95) 
(A 96) 
(A 97) 
(A 98) 
(A 99) 
(A 100) 
(A 101) 
(A 102) 
(A 103) 
(A 104) 
(A 105) 
(A 106) 
(A 107) 
(A 108) 
(A 109) 
(A 110) 
(A 111) 
(A 112) 
(A 113) 
(A 114) 
(A 115) 
(A 116) 
(A 117) 
(A 118) 
(A 119) 
(A 120) 
(A 121) 
(A 122) 
(A 123) 
(A 124) 
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IF (GAHMA.LT. 2. 0005) GO TO 41 


(A 

125) 


GAMMA=2 • 0 


(A 

126) 


GO TO 812 


(A 

127) 

41 

CONTINUE 


(A 

128) 


IF ( ITAIL.EQ.O) GO TO 1 


(A 

129) 


SUML=0. 031526 


(A 

130) 


IF (LI.EQ.l) SUML=0 .025 


(A 

131) 


IF (S IGMAW.EQ.O. .AND. S I GMAT. EQ .0 . ) 

GO TO 450 

(A 

132) 


IF ( S IGMAW.EQ.O. . AND. S I GMAT . NE .0 . ) 

GO TO 455 

(A 

133) 


IF (SIGMAW.NE.O. .AND. SI GMAT. EQ.O . ) 

GO TO 460 

( A 

134) 


M 7=4 


(A 

135) 


N7 = 10 


(A 

136) 


CONST 1 = 1 • 0 


(A 

137) 


IF (ETA1.NE.1.) GO TO 412 


(A 

138) 


M7=2 


(A 

139) 


CONST 1=2 .0 


(A 

140) 


GO TO 412 


(A 

141) 

450 

M7=N7=1 


(A 

142) 


XLOAD(1)=1.0 


(A 

143) 


SUML=0. 02 5 


(A 

144) 


CONST 1=40 . 0 


(A 

145) 


GO TO 412 


(A 

146) 

455 

M 7 = 4 


(A 

147) 


N7 = l 


(A 

148) 


XLOAD ( 1 ) = 1 .0 


(A 

149) 


SUML = 0 .02 5 


( A 

150) 


CONST1=10.0 


(A 

151 ) 


IF (ETAi.NE.l.) GO TO 412 


(A 

152) 


H7=2 


(A 

153) 


CONST 1=20. 0 


(A 

154) 


GO TO 412 


(A 

155) 

460 

M7 = l 


(A 

156) 


N7=10 


(A 

157) 


CONST 1 = 4 . 0 


(A 

158) 


IF (ETAi.NE.l.) GO TO 412 


(A 

159) 


N7=5 


(A 

160) 


CONST 1 = 8.0 


(A 

161) 

412 

DO 420 I TL= 1 , I T A I L 


(A 

162) 


WRITE (6,901) SIGMA, I AL PHA , AL AM , AZ ET A , E TA 1 , AALP 

( A 

163) 


WRITE (6,210) 


( A 

164) 


WRITE (6,211) 


(A 

165) 


WRITE (6,212) 


(A 

1661 


WRITE (6,213) 


(A 

167) 


WRITE (6,214) 


(A 

168) 


WRITE (6,215) 


(A 

169) 


WRITE (6,216) 


( A 

170) 


WRITE (6,217) 


(A 

171) 


WRITE (6,218) 


( A 

172) 


DO 91 K= 1,8 


(A 

173) 


GAMMA=STOR E ( K ) 


(A 

174) 

399 

AC=0. 5*SIGMA*TAN(LAMBDA) *GAMMA 


(A 

175) 


IF (LI.NE.l) AC=0.848*AC 


(A 

176) 


TL = ( AT L* FLOAT ( I TL )* GAMMA/ 2 .0 ) +AC 


(A 

177) 


TH=ATH*FLOAT ( ITL ) *G AMMA/ 2 . 0 


(A 

178) 


ZETAl=AZETA/( 1.0 +AC*AZETA*SIN( ALPHA) ) 

(A 

179) 


ZETAl=2.0/(2.0+GAMMA*( 1.0-ZETA1) ) 


( A 

180) 


DO 400 L 2=1 , 28 


( A 

181) 


XDELT A ( L2 ) =0 . 


(A 

182) 

400 

DELTA ( L2 ) =0 . 


(A 

183) 


DO 401 M 1 = 1 , M7 


(A 

184) 


DO 402 Nl= 1 , N7 


(A 

185) 
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XSTAR= ( 2. AFLOAT ( Nil -11. 1/10,0 
YSTAR=( 11.-2. *FL0AT(N1) ) /10.0 
ZSTAR=(5.-2.*FLOAT( Ml) )/4.0 
ETA=ETA1+XSTAR* SIGMA 

ZETA=ZETAl/( l.-ABSI YSTAR ) *S I GMA*GAMMA*ZE TA1*T AN( L AMBDA )* 

1 SIN ( ALPHA) ) 

XOVERH=TL*COS( ALPHA >+TH*S INC AL PHA >— SI GMA *GAMMA*T AN ( LAMBDA ) * 
1 COS( ALPHA) *ABS (YSTAR ) 

YOVERH=ZSTAR*SIGMAT*GAMMA-YSTAR*SIGMA*GAMMA 

ZOVERH=TH*COS( ALPHA )-TL*SIN< ALPHA ) +S IGMA*GAMMA*T AN ( L AMBDA ) * 
1 SZN( ALPHA) *ABS( YSTAR) 

CALL DLTAS (C(K)) 

********************* SEE APPENDIX G FOR SUBROUTINE DLTAS 
DO 405 L 1= 1 , 2 8 

405 XDELTAUl ) =XDEL T A ( L 1 ) + 1 DEL TA { L I ) *XLOAD( N 1 ) ) 

402 CONTINUE 
401 CONTINUE 

DO 407 L 3= 1 * 2 8 

407 DELTAIL3) -XDELT A(L3)*SUML*CQNST1*( GAMMA/ 2.0) 

TL=ATL* FLOAT ( ITL ) 

TH=ATH*FLOAT ( ITL ) 

WRITE (6,148) C(K) ,GAMM A , SI GMAT , TL , TH 
WRITE (6,150 (DELTA(I) ,1=1,25,4) 

WRITE (6,151) ( DELTA! I) ,1=2, 26, A) 

WRITE (6,152) ( DELTA! I) ,1=3*27,4) 

WRITE (6,153) (OELTA( I) ,1=4,28,41 
DO 414 L4=l ,28 
XDELTA(L4)=0-0 
IF ( GAMMA. LT. 2. 0005) 

GAMMA=2 • 0 
GO TO 399 
CONTINUE 
CONTINUE 
GO TO 1 

FORMAT (5X*IMPINGEMENT 
FORMAT (5X*IMPINGEMENT 


(A 

186) 

(A 

187) 

(A 

188) 

(A 

189) 

(A 

190) 

( A 

191) 

(A 

192) 

( A 

193) 

(A 

194) 

(A 

195) 

( A 

196) 

(A 

197) 


********************** 


414 


91 

420 


GO TO 91 


119 

120 


150 

151 

152 

153 
210 
211 


IS INFINITE*) 
(SPANS) USING CHI 


(EFFECTIVE) IS* 


FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 


DISTANCE 
DISTANCE 

1 FI 0 .2 , 5X*U SING CHI (MOMENTUM) IS*F10.2) 

148 FORMAT <//5H CH I =F 5 .1 , 5X6HGAMMA=F7 . 3 , 5X* I NT ERFERENCE AT TAIL*5X 
1 *S I GMA (TAIL)=*F6*3, 5X*TL /H=*F7 . 3 , 5X*TH/H=*F7 . 3/ ) 

149 FORMAT (//5H CH I =F 5 . 1 , 5X6HGAMMA= F7. 3 , 5X* INT ERF ER ENC E AT WING*/) 
(3X5H{ W,L ) 7( F17.4) ) 

(3X5H(U,L )7(F17.4) ) 

( 3 X 5H ( W, D ) 7 ( F17.4) ) 

(3X5HCU, D)7(F17«4)) 

( 1X131(1 H- ) ) 

(1X1HI11X 1HI31X61HC0RRECTICN FACTORS FOR CORRECTING FROM A 
1 WIND TUNNEL WHICH IS25X1HI) 

212 FORMAT ( 1 X 1H I 1 1 X 1H 1 1 1 7 ( 1 H- ) 1 H I ) 

213 FORMAT ( 1 X 1H I 11 X 1H l 16X 1H I 5X6HCL0S ED5X 1H I 16X1 H I 2X1 2HCL0S ED FLOOR2X1 
1 H I 6X4 HOP EN 6 X 1HI 1 6X1 HI 5X6HCL0SED4X1HI ) 

214 FORMAT ( 1 X 1H I 3X 5 HDELTA3X 1HI 5 X6HCLCSED5X1 H I 4X9H0N B0TT0M3X1H I 6X4H0P 
1EN6X1HI6X4HONLY6X1HI5X5HFLOOR6X1HI 5X6HCL 0SE05X1HI 3X9 HON B0TT0M3X 1H 
21 ) 

215 FORMAT ( 1 X1H I 11 X 1H I 1 6X1H 1 6X4H0NL Y6X 1H I 16X 1 8H I ( GROUND EFFECT) I6X4H 
1GNLY6X1HI 16X1HI6X4HCNL Y5X1HI ) 

216 FORMAT ( 1X1HI 11X 1HI 84< 1H-) IHI32C 1H- )1HI ) 

217 ■ FORMAT i IX 1HI 11X 1H I 36X11HTO FREE AI R37X1 HI 8X16HTO GROUND EFFECT8X1 
1 H I ) 

218 FORMAT 

900 FORMAT 

901 FORMAT 


(1X131(1 H- ) ) 

( II ,F9.3, 4F10.3/I1,F9.3,2F10.3) 

( 1H1///40X*AVERAGE INTERFERENCE OVER A SWEPT 


OF FINITE 
CLOSED-ON 


999 


WING 

1 S PAN*//30X*GAMM A CCMPUTED TO YIELD ZERO DEL T A ( W , L ) IN A 
2-BOTTOM-ONLY TUN NEL*///2 9X*S IGMA = *F6 . 3 , 16X , A8* LOAD ING + 15X, 
3*LAMBDA =* F7 • 3/ / 29 X *Z ETA =*F6 . 3 , 18X*ETA -*F7. 3, 17X* ALPHA =* 
4F7.3///35X*ALL VALUES ARE REFERENCED TO THE HIGH-SPEED ( GAMM A= 2.0) 

5 SECT I0N*/39X*ALL DIMENSIONS ARE MEASURED FROM THE AERODYNAMIC CEN 

6 TER*//) 

STOP 

END 


( A 

198) 

( A 

199) 

(A 

200) 

( A 

201 > 

< A 

202 ) 

( A 

203) 

(A 

204) 

( A 

205) 

(A 

206) 

< A 

207) 

(A 

208) 

(A 

209) 

(A 

210) 

(A 

211) 

( A 

212) 

(A 

213) 

(A 

214) 

(A 

215) 

(A 

216) 

( A 

217) 

(A 

218) 

(A 

219) 

( A 

220 ) 

(A 

221) 

(A 

222) 

(A 

223) 

( A 

224) 

( A 

225 ) 

( A 

226) 

(A 

227) 

( A 

228) 

( A 

229) 

(A 

230) 

( A 

231) 

(A 

232) 

(A 

233) 

(A 

234) 

( A 

235) 

(A 

236) 

(A 

237) 

(A 

238) 

( A 

239) 

( A 

240) 

(A 

241 ) 

( A 

242) 

(A 

243 > 

( A 

244) 

( A 

245) 

( A 

246) 

( A 

247) 

(A 

248) 

( A 

249) 

( A 

250) 

(A 

251) 

(A 

252) 

(A 

253) 
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FORTRAN PROGRAM FOR COMPUTING INTERFERENCE OVER ROTORS 
IN A VARIABLE WI OT H— HE IGHT— RAT 1 0 CLOSED-ON-BOTTOM-ONLY TUNNEL 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN, VERSION 2.1, TO RUN ON CDC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


NOTE THAT ALL INPUT AND OUTPUT VARIABLES ARE REFERENCED TO THE GAMMA*2.0 
TEST SECTION. INPUT IS AT ADDRESS 1 IN FORMAT 900 I TWO CARDS PER CASE). THE 
REQUIRED INPUT VARIABLES ARE 


LI DISK-LOAD DISTRIBUTION INDICATOR, LI»1 FOR UNIFORM LOAOING, 

LI =2 FOR TRIANGULAR LOAOING 

ZETA1 SEMIHEIGHT OF GAMMA-2 SECTION OIVIOED BY HEIGHT OF ROTOR 

CENTER ABOVE FLOOR IN GAMMA-2 SECTION 

ETA1 DISTANCE FROM ROTOR CENTER TO RIGHT-HAND WALL OIVIDEO 

BY TUNNEL SEMIMIOTH 

SIGMA RATIO OF ROTOR DIAMETER TO TUNNEL WIDTH 

ALPHA ANGLE OF ATTACK OF ROTOR TIP-PATH PLANE, DEG 

I TAIL NUMBER (0 TO 9) OF TAIL POSITIONS, PROPORTIONATELY INCREASED 

IN TAIL LENGTH AND HEIGHT, AT WHICH INTERFERENCE IS 
REQUIREO. USE SMALLEST LENGTH AND HEIGHT AS INPUT. 

S IGNAT RATIO OF TAIL SPAN TO TUNNEL WIOTH 

TL DISTANCE OF TAIL BEHIND ROTOR CENTER, NONDIMENSIONALIZED WITH 

RESPECT TO ROTOR RADIUS 

TH HEIGHT OF TAIL ABOVE ROTOR CENTER, NONDIMENSIONALIZED WITH 

RESPECT TO ROTOR RADIUS 

ALPHAB ANGLE OF ATTACK OF BODY CARRYING TAIL, DEG 


THIS PROGRAM REQUIRES SUBROUTINE DLTAS (SEE APPENDIX G) AND SUBROUTINE 
VARLET (SEE APPENDIX HI. 


SINCE THE TAIL LOCATIONS ARE NCNOIMENS I ONAL I ZED WITH RESPECT TO ROTOR 
RADIUS, THIS PROGRAM CAN NOT ACCOMMODATE, ANO WILL REJECT, TAIL INTERFERENCE 
CALCULATIONS FOR CASES INVOLVING SIGMA * 0, HOWEVER, THE INTERFERENCE AT THE 
CENTER OF LIFT WILL BE CALCULATED FOR SUCH CASES. TAIL INTERFERENCE FOR SUCH 
CASES CAN BE OBTAINED FROM THE CORRESPONDING PROGRAM FOR WINGS, SINCE, WHEN 
SIGMA IS ZERO, THE REPRESENTATION OF WINGS AND ROTORS IS IDENTICAL. 
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PROGRAM W I NDTUN ( INPUT, OUT PUT , TAP E5= INPUT , TAP E6=OUTPU T ) 

COMMON ZETA, ETA, GAMMA , XOVERH , YOVERH , ZOVE R H , DELTA ( 28 ) 

DIMENSION XDELTA(28),PSI(20) ,XLQAD(2G) , RUN I F ( 20 ) , RTR I A t 20 ) , 

1 C(8) , STORE! 8 ) 

DATA ( RUN I F ( I ) , 1=1 , 20)/4*0.2981, 8*0. 6255, 8*0-8921/ 

DATA ( RTR I A ( I) , I =1 , 20 ) /4*0 .4 38 6 , 8*0. 7296, 8*0 . 9262 / 

DATA <C(I) , 1=1,8 >/2G.,30.,40., 50. ,60. ,70., 80., 90./ 

R AD=0 .0174532925199 
PI=3. 14159265358979 
PS I ( 1 ) = ( P I /4 • ) 

PS I ( 2 ) = 3 • *PS I ( 1 ) 

PSl(3)=5.*PSI(l) 

PS I ( 4 ) = 7 • *PS I ( 1 ) 

PSI ( 5 )-PS I 1 13 ) = ( PI/8. » 

PSl(6)=PSI(14)=3.*PSI(5) 

PSI<7)=P$I(15)=5.*PSI(5) 

PS I ( 8 ) = PS I ( 1 6 ) = 7 • *P S I t 5 ) 

P S I < 9 ) = P S I ( 1 7 > = 9 . *P S I ( 5 ) 

P$l( 10) = PSI( 18) = 11.*PSI( 5) 

PSI(11)=PSI (19)=13.*PSI(5) 

PSI( 12)=PSI( 20)= 15. *PSI(5 ) 

1 READ (5,900) L I ,ZE TA1 ,ETA1, SI GMA , ALPHA, I TA I L, SI GMAT ,TL , TH, ALPHAS 
IF (EOF, 5) 999,47 
47 IF (LI.EQ.l) GO TO 804 
I ALPHAS GHTR I ANGULAR 
DO 808 M2=l,20 

808 XL0AD(M2)=RTRIA( M2) 

GO TO 160 

804 I ALPH A= 10 H UNIFORM 
DO 809 M 2= 1 , 20 

809 XL0AD(M2)=RUNIF( M2) 

160 WRITE (6,901) I ALPHA , SI GMA , ALPH A, Z ETA1 , SIGMAT , ALPHAS, ETA1 
WRITE (6,210) 

WRITE (6,211) 

WRITE (6,212) 

WRITE (6,213) 

WRITE (6,214) 

WRITE (6,215) 

WRITE (6,216) 

WRITE (6,217) 

WRITE (6,218) 

SUML=0 .0025 
CONST 1 = 1 . 

AZETA=ZETA1 
AT L=T L 
ATH= TH 
AL PR- ALPHA 
ALPB= ALPHAS 
ALPHA= ALPR *R AD 
ALPHAB=ALPB*RAD 
DO 41 K= 1 , 8 
DO 803 L 1= 1 , 2 8 
DE LT A ( L 1 ) =0 . 

803 XDELT A( LI ) =0 . 

I NK = 1 

GAMMA=GM Al=2 • 0 
M 7 = N7 = 20 

IF (SIGMA. NE.O. ) GO TO 811 


(B 1) 
(B 2) 
(B 3) 
(B 4) 
(B 5) 
(B 6) 
(B 7) 
(B 8) 
(B 9) 
(B 10) 
(B 11) 
(8 12 ) 
(B 13) 
(B 14) 
(B 15) 
(B 16) 
(B 17) 
(B 18) 
(B 19) 
(B 20) 
(8 21 ) 
(B 22) 
(B 23) 
(B 24) 
(B 25) 
(B 26) 
(B 27) 
(B 28) 
(B 29) 
(B 30) 
(B 31) 
(B 32) 
(B 33) 
(B 34) 
(B 35) 
(B 36) 
(B 37) 
(B 38) 
(B 39) 
(B 40) 
(B 41) 
(B 42) 

( B 43) 
(B 44) 
(B 45) 
(B 46) 
(B 47) 
(8 48) 
(B 49) 
(B 50) 
(B 51) 
(B 52) 
(B 53) 
(B 54) 
(B 55) 
(8 56) 
(B 57) 
(B 58) 
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N7=M7= 1 


( B 

59) 


CONST 1— 4G0 


(B 

60) 


GG TO 812 


(B 

61 ) 

811 

IF (ETA1.NE.1.) GO TO 812 


(B 

62) 


CONST 1 = 2 • 


(B 

63) 

812 

DO 801 M 1 = 1 , M7 


( B 

64) 


DO 802 N1 = 1 » N7 


(B 

65) 

836 

IF (ETA1.NE.1.) GO TO 840 


(B 

66) 


IF (PSI C N 1 ).GT.P I) GO TO 802 


(B 

67) 

840 

ETA = ETAl-( XLO AD ( N1 )*SIGMA *S I N ( P S I < N1 ) ) ) 


(B 

68) 


ZETAl=2.0/(2.0+GAMMA*( 1.0-AZETA) ) 


(8 

69) 


ZETA=1 • / ( ( 1./ZETA1)-CXL0AD(N1 )*S IGMA *S IN ( ALPHA ) *COS < PS I ( N 1 ) ) *GAMM 

(B 

70) 

1A) ) 


(B 

71) 


XOVERH-S IGMA *GAMMA*COS< ALPHA)*! XLO AD ( Ml )*COS(PSI (Ml ) >-XLOAD ( N1 ) *C 

(B 

72) 

10S(PSI(N1 ) ) ) 


( B 

73) 


YOVERH= SIGMA *GAMM A* ( XLO AD ( M 1) *S IN (PS I (Ml) )- XLOAD( N1 ) *S IN ( PS I ( N1 ) ) 

(B 

74) 


1 ) 


(B 

75) 


ZOVERH— SIGMA *G AMM A*S IN ( ALP HA ) * ( XLO AD( M 1 ) *COS ( P S I ( M 1 ) )-XLOAD(Nl )* 

(B 

76) 


1C0S ( P S I ( N 1 ) ) ) 


(B 

77) 


IF ( INK.EQ.7 ) GO TO 850 


(B 

78) 


CALL VARLET C CC K ) 1 


(B 

79) 

********************* SEE APPENDIX h for subroutine 

VARLET ********************* 


GO TO 832 


(B 

80) 

850 

CALL DLTAS (C(K)) 


(8 

81) 

********************* SEE APPENDIX G FOR SUBROUTINE 

DLTAS ********************** 

832 

DO 805 L 1 = 1 » 2 8 


(8 

82) 

805 

XDELTA(L1 )=XDELTA( L1)+(DELTA( LI ) ) 


( ft 

83 ) 

802 

CONTINUE 


( B 

84) 

801 

CONTINUE 


( B 

85 ) 


IF ( INK.EQ.7 ) GO TO 835 


( B 

86) 


DELTA! 5) = XDELTA( 5 ) *SUML*CONS Tl*( GAMMA/2.0) 


( B 

87) 


IF { ABS ( DELTA(5) ) .LT. 0.00005) GO TO 833 


(B 

88) 


DEL2=DELTA(5> 


(B 

89) 


IF ( INK.NE.l ) GO TO 831 


(B 

90) 


DEL 1=DELT A ( 5 ) 


(B 

91 ) 


GAMMA- GMA2-1.0 


(B 

92) 


INK-2 


(B 

93) 


XDELTA(5)=0.0 


(B 

94) 


GO TO 812 


(8 

95) 

831 

IF (ABS(DELTA(5) ).LT. 0.00005) GC TO 833 


(B 

96) 


GAMMA-GAMM A- ( GMA 2- GMA1 ) * DEL 2 / ( DEL2-DEL 1 ) 


( B 

97) 


IF (GAMMA.LE.O. 1 i GAMMA -0.1 


(B 

98) 


DEL1-DEL2 


(B 

99) 


GMA1-GMA2 


( B 

100) 


INK-3 


( B 

101 ) 


GM A2-GAMMA 


( B 

102 ) 


IF (GMA2.EQ.GMA1 ) GO TO 837 


( B 

103) 


XDELTA(5)=0.0 


( B 

104) 


GO TO 812 


( 8 

105) 

837 

STORE ( K ) -0 .01 


( B 

106) 


GO TO 41 


( B 

107) 

833 

I NK = 7 


( B 

108) 


XDEL T A ( 5 ) = 0 • 0 


( B 

109) 


STORE (K)-GAMMA 


( B 

110) 


GO TO 812 


( B 

111 ) 

835 

DO 807 L3 = l f 28 


(B 

112 ) 

807 

DELTA(L3 )=XDELTA( L3 )*SUML*C0NST1* (GAMMA/2. ) 


( B 

113) 


WRITE (6,149) C(K), GAMMA 


(B 

114) 


WRITE (6,150) (DELTA! I), 1 = 1,25,4) 


( B 

115) 


WRITE (6,151) (DELTA! I), 1 = 2,26,4) 


( B 

116) 


WRITE (6,152) (DELTA! I), 1 = 3,27,4) 


( B 

117) 
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WRITE (6*153) <DELTA( I)* 1=4*28, 4) 

IF ( S IGMA • EQ.O. ) GO TO 822 
IF <C(K).EQ.90.> GO TO 820 
XS1=TAN(C( K)*RAD )/( 2 .0 *S I GMA *GAMM A*ZETA l ) 

CHIM=2.0*C(K)-90. 

XS2=TAN(CHIM*RAD ) / ( 2. 0 *S I GMA*GAMMA*ZETA1 ) 

WRITE (6*120) X SI , XS2 
GO TO 822 
820 WRITE (6,119) 

822 DO 814 L4= 1,28 
814 XDELT A ( L4 ) =0 . 

IF (GAMMA.LT. 2. 0G05) GO TO 41 
GAMMA=2.Q 
GO TO 812 
41 CONTINUE 

IF ( ITAIL.EQ.O) GO TO 1 
IF (SIGMA. NE. 0. ) GO TO 400 

WRITE (6*901 ) IALPHA,SIGMA,ALPR*AZETA»SI GMAT *ALPB*ETA1 
WRITE (6,904) 

GO TO 1 

400 SUML = 0. 01 2 5 

IF (SIGMAT.NE.O. ) GO TO 411 
M7=l 

CONST 1=4 .0 
GO TO 412 

411 IF (ETA1.NE.1.) GO TO 413 
M7=2 

CONST 1=2 .0 
GO TO 412 
413 M7 =4 

CONST 1 = 1 .0 

412 00 420 ITL=1,ITAIL 
TL=ATL*FLOAT ( ITL ) 

TH=ATH*FLOAT ( ITL ) 

WRITE (6,901) I ALPHA,SIGMA*ALPR,AZETA,SIGMAT,ALPB,ETA1 
WRITE (6,210) 

WRITE (6,211) 

WRITE (6,212) 

WRITE (6,213) 

WRITE (6,214) 

WRITE (6,215) 

WRITE (6,216) 

WRITE (6,217) 

WRITE (6,218) 

DO 91 K= 1 , 8 
GAMMA=STORE(K) 

399 ZETAl=2.0/(2.0+GAMMA*(1.0-AZETA) ) 

DO 401 Ml =1 , M7 
DO 402 N 1=1 , 20 

ETA=ETA1-(XL0AD( N1)*SIGM A*S I N ( PS I ( N1 ) ) ) 

ZETA=l./< ( 1./ZETA1)-XL0A0(N1)*SIGMA*GAMMA*SIN( ALPHA) *COS(PSI (N1 ) ) ) 
XOVERH=SIGMA*GAMMA*( (TL*COS( ALPHAS) )+(TH*S IN (ALPHAS ) )-(XLOAD(Nl) 
l *C0S ( ALPHA) *C0S(PSI(N1)) ) ) 

XM1=FL0AT ( Ml ) 

YOVERH=SIGMA*GAMMA*(-< (2.*XMl-5. )/4. ) *( S I GMAT/ S I GMA ) - ( XLOAD( N1 ) 
l *SIN ( PSI (N1 ) ))) 

ZOVERH=SIGMA*GAMMA*( ( TH* COS ( ALPHAB > )-( TL *S IN ( ALPHAS ) )+(XLOAD(Nl) 

1 *SIN( ALPHA)*C0S(P$I(N1 >) ) ) 

CALL DLTAS (C(K)) 


(B 118) 
(B 119) 
(B 120) 
(B 121) 
(B 122) 
(B 123) 
(B 124) 
(B 125) 
( B 126) 
(B 127) 
(B 128) 
(B 129) 
(B 130) 
(B 131) 
(B 132) 
(B 133) 
(B 134) 
(B 135) 
(B 136) 
(B 137) 
(B 138) 
(B 139) 
(B 140) 
(B 141) 
(B 142) 
(B 143) 
(B 144) 
(B 145) 
(B 146) 
(B 147) 
(B 148) 
(B 149) 
(B 150) 
(B 151) 
(B 152) 
(B 153) 
(B 154) 
(B 155) 
(B 156) 
(B 157) 
(B 158) 
(B 159) 
(B 160) 
(B 161) 
(B 162) 
(B 163) 
(B 164) 
(B 165) 
(B 166) 
(B 167) 
(B 168) 
(B 169) 
(B 170) 
(B 171) 
(B 172) 
(B 173) 
(B 174) 
(B 175) 

( B 176) 


********************* SEE APPENDIX G FOR SUBROUTINE DLTAS ********************** 


DO 405 L 1= 1 , 28 


(B 177) 
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405 XDELTAILl )=XDELTA(L1)+DELTA(L1) ( B 

402 CONTINUE (B 

401 CONTINUE (B 

DO 407 L 3= 1,28 (B 

407 DELTA(L3) = XDELTA (L3 ) *SUML*CONST 1 * ( GAMMA/ 2.0) ( B 

WRITE (6,148) C ( K I , GAMMA , TL , TH (B 

WRITE (6,150) (DELTA( I ) , 1=1,25,4) (B 

WRITE (6,151) (DELTA(I), 1=2,26,4) (B 

WRITE (6,152) (DELTA(I), 1=3,27,4) (B 

WRITE (6,153) (DELTA(I ) , 1=4,28,4) (B 

DO 414 L4= 1 , 28 ( B 

414 XDELTA(L4)=0.0 (B 

IF (GAMMA. LT. 2. 0005) GO TO 91 (B 

GAMMA=2. 0 (B 

GO TO 399 (B 

91 CONTINUE (B 

420 CONTINUE (B 

GO TO 1 (B 

119 FORMAT ( 5X*IMPI NGEMENT DISTANCE IS INFINITE*) (B 

120 FORMAT ( 5 X* I MP I NGE MENT DISTANCE (SPANS) USING CHI (EFFECTIVE) IS* (B 

1F10.2,5X*USING CHI (MOMENTUM) IS*F10.2) (B 

148 FORMAT (//5X*CHI =*F7 . 3, 5X+GAMMA =*F 7 . 3, 5 X* I N T ERF ERENC E AT TAIL*5X (B 

1*TL/R =*F6.3,5X*TH/R =*F6.3/) (B 

149 FORMAT (//5H CH l =F 5 . 1 , 5X 6HGA M MA = F7 . 3 , 5 X* I NT ER F ER ENC E AT ROTOR*/) (B 

150 FORMAT ( 3 X 5H ( W, L ) 7 ( F 17 .4 ) ) (B 

151 FORMAT (3X5H(U,L )7( F17.4) ) (B 

152 FORMAT ( 3X5H ( W, D ) 7 ( F17 .4 ) ) (B 

153 FORMAT ( 3 X 5H ( U, D ) 7 ( FI 7 .4 ) ) (B 

210 FORMAT (1X13K1H-)) (B 

211 FORMAT (1X1HI 11X 1HI31X61HCQRRECTICN FACTORS FOR CORRECTING FROM A (B 

1WIND TUNNEL WHICH IS25X1HI) (B 

212 FORMAT ( 1 X 1H 1 11 X 1H I 1 1 7 ( 1 H- ) 1 H I ) (B 

213 FORMAT ( 1 X 1H 1 11 X 1H I 1 6X 1H I 5X6HCL0 SED5X 1H I 1 6 X 1H I 2X 1 2HC LOSED FL OOR2 XI (B 

1 H I 6X4H0PEN6 X 1HI 1 6X1 HI 5X6HCL0SED4X 1H I ) ( B 

214 .FORMAT ( 1 X1HI 3X 5HDELTA3X 1HI5X6HC LCSED5X1 H I4X9H0N BOTTOM 3X1HI6X4 HOP (B 

1EN6X1HI6 X4H0NLV6 X1H I5X5HFL00R6X1HI 5X6 HCL 0SED5X1HI 3X9H0N BOTT QM3 X 1H ( B 
21 ) (B 

215 FORMAT ( 1 X 1 H 1 11 X 1H 1 16X 1H I6X4H0NLY6X 1 HI 16 XI 8H I ( GROUND EFFECT) I6X4H (B 

1GNLY6X1HU6X1HI6X4HCNLY5X1HI ) (B 

216 FORMAT ( 1 X 1H 1 11 X 1H I 84 ( 1H- ) 1H 1 3 2( 1H- ) 1H I ) (B 

217 FORMAT ( 1 X 1 H 1 11 X 1 H I 3 6X 1 1 HTO FREE A l R 37 X 1 H I 8X 1 6HT 0 GROUND EFFECT8X1 (B 

1 H I ) (B 

218 FORMAT (1X13K1H-)) (B 

900 FORMAT ( 1 1 ,F 9 . 3 , 3F 1 C . 3/ 1 1 ,F9 . 3 ,3F 10 . 3 ) (B 

901 FORMAT ( 1 H 1/ / // 4 2X * AV E RA GE INTERFERENCE OVER A ROTOR OF FINITE SPA (B 

1N*//30X*GAMMA COMPUTED TO YIELD ZERO DELTA ( W,L ) IN A CLOS ED- ON-BOT (B 
2 TOM-ONLY T UNNEL * // / 55 X , A 10* DISK L OA D I NG */ / 30 X*S I GMA ( ROTOR ) =* (B 

3F6.3, 8X*ALPHA(R0T0R ) =*F 7*3, 10X*ZETA =*F6.3//30X*SIGMA( TAIL) =* (B 

4F6.3, 8X*ALPHA(TA IL ) =*F 7*3, 10X*ETA =*F 6 • 3/ / 3 5X* ALL VALUES ARE RE (B 

5F ERENC ED TO THE HIGH-SPEED (GAMMA=2.0) SECTION*//) (B 

904 FORMAT {///40X*TAIL IS OMITTED - EQUATIONS ARE NOT VALID WHEN SIGM (B 
1 A ( ROT OR ) = 0.*///) (B 

999 STOP (B 

END (B 


178) 

179) 

180) 
181 ) 
182) 

183) 

184) 

185) 

186) 

187) 

188) 

189) 

190) 

191) 

192) 

193) 

194) 

195) 

196) 

197) 

198) 

199) 

200 ) 
201 ) 
202 ) 

203) 

204) 

205) 

206) 

207) 

208) 

209) 

210 ) 
211 ) 
212 ) 

213) 

214) 

215) 

216) 

217) 

218) 

219) 

220 ) 
221 ) 
222 ) 

223) 

224) 

225) 

226) 

227) 

228) 
229) 
2 30) 
231) 
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FORTRAN PROGRAM FOR CCMPUT ING INTERFERENCE OVER WINGS 
IN A VARIABLE MODEL-HEIGHT CLOSED- ON-BOTTOM-ONLY TUNNEL 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN , VERSION 2,1, TO RUN ON CDC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROX I MAT EL Y 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DCUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


INPUT IS AT ADDRESS 1 IN FORMAT 9CC (TWO CARDS PER CASE). THE REQUIRED 
INPUT VARIABLES ARE 


LI SPAN LOAD- 01 STR I BUT ICN INDICATOR, LI = I FOR UNIFORM LOADING, 

L I =2 FOR ELLIPTIC LOADING 

ETA1 DISTANCE FROM AERODYNAMIC CENTER TO RIGHT-HAND WALL OIVIDEO 

BY TUNNEL SEMIWIOTH 

GAMMA WIDTH-HEIGHT RATIO OF WIND TUNNEL 

SIGMA RATIO OF WING SPAN TO TUNNEL WIDTH 

LAMBDA WING SWEEP ANGLE, DEG 

ALPHA ANGLE OF ATTACK OF WING, DEG 

ITAIL NUMBER (0 TO 9) OF TAIL POSITIONS, PROPORTIONATELY INCREASED 

IN TAIL LENGTH AND HEIGHT, AT WHICH INTERFERENCE IS 
REQUIRED. USE SMALLEST LENGTH AND HEIGHT AS INPUT. 

SIGMAT RATIO OF TAIL SPAN TO TUNNEL WIDTH 

TL DISTANCE OF TAIL BEHIND AERODYNAMIC CENTER, NON D I MENS IONAL I ZE D 

WITH RESPECT TO TUNNEL SEMIHEIGHT 

TH HEIGHT OF TAIL ABOVE AERODYNAMIC CENTER, NO ND I ME NS I ONAL I Z E D 

WITH RESPECT TO TUNNEL SEMIHEIGHT 


THIS PROGRAM REQUIRES SUBROUTINE DLTAS (SEE APPENDIX G) AND SUBROUTINE 
VARLET (SEE APPENDIX H). 


PROGRAM WINDTUN( INPUT, OUTPUT , T AP E5= I NPUT , TAP E6=0UTPU T ) (C 1) 
COMMON ZETA, ETA, GAMMA, XOVERH,YOVERH,ZOV£RH,DELTA(28) (C 2) 
DIMENSION XLE(10 ) ,XLOAD( 10) , XDELTA(28) ,C( 8) ,ST0RE(8 > (C 3) 
REAL LAMBDA (C 4) 
DATA t C ( I ) ,1-1,8 )/20. ,30. ,40. ,50. ,60., 70. ,80. ,90./ (C 5) 
XLE(1 )=XLE(10)=0. 43579 (C 6) 
XLE(2)~XLE(9 )~0 • 71422 (C 7) 


38 
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l 

47 

808 

804 

809 

160 


803 


811 

813 

812 


XLE(3)=XLE(8)=0. 86603 
XLE(4)=XLE (7 ) = 0 . 95394 
XLE(5)=XLE(6) =0.99499 

READ ( 5 » 900 ) LI , ET A1 , GAMMA , S I GMA , LAMBDA , ALPHA, ITAILtSI GMAT ,TL,TH 

IF (EOF, 51 999,47 

IF (LI.EQ.l) GO TO 804 

IALPHA=8HELLIPTIC 

SUML=0. 0126104 

DO 808 M2= 1,10 

XLOAD(M2)=XLE(M2 ) 

GO TO 160 
SUML =0 .01 
I ALPHA=8HUNIFORM 

00 809 M2= 1 , 10 
XLOAD ( M2 ) = 1 • 0 

WRITE (6,901 ) SIGMA, I ALPHA, LAMBDA, GAMMA, ETA1 , ALPHA 
WRITE (6,210) 

WRITE (6,211) 

WRITE (6,212) 

WRITE (6,213) 

WRITE (6,214) 

WRITE (6,215) 

WRITE (6,216) 

WRITE (6,217) 

WRITE (6,218) 

AA LP = ALPHA 

AL AM= L AMBDA 

AT L = T L 

ATH=TH 

CONST 1=1.0 

RAD=0. 0174532925199 

ALPHA= ALPHA* RAD 

LAMBDA=LAMBDA*RAD 

AC= 0.5*SIGMA*GAMMA*TAN( LAMBDA) 

IF (LI.NE.l) AC = 0 • 848*AC 
DO 41 K= 1,8 
DO 803 L 1 = 1 , 28 
DELTA(L1 )=0. 

XDELT A ( L 1 ) =0 • 

1 NK = 1 

ZETA1=ZET1=1.0 

IF (SIGMA. NE. 0. ) GO TO 811 

M7=N7 = 1 

XLOAD ( l ) = 1 . 0 

SUML= 1 • 

GO TO 812 

IF (ETA1.NE.1.) GO TO 813 
M7 = 5 
N7 = l 0 
CONST 1 = 2 . 

GO TO 812 

M7=N7=10 

DO 801 M1=1,M7 

DO 802 N 1 =1 , N7 

X$TAR=(11. —2 • *FL OAT ( M 1 ) ) / 10 • 

YSTAR=(2.*FL0ATC ND-l 1.1/10. 

ZSTAR=(11.-2.*FL0AT (Nil ) /10. 

E T A=ET A1 +YSTAR* SIGMA 

ZETA=ZETA1/ ( 1 .-A BS ( YSTAR ) *S I GMA* GAMMA *ZE TA1*TAN( LAMBOA)*SIN( ALPHA) 
1 ) 

XO VERH= SI GMA*GAM MA*T AN( L AMBDA )*C OS ( ALPHA )*(ABS(XSTAR )-ABS(ZSTAR) ) 


(C 8) 
(C 9) 
(C 10) 
(C 11) 
(C 12) 
(C 13) 
(C 14) 
(C 15) 
(C 16) 
(C 17) 
(C 18) 
(C 19) 
(C 20) 
(C 21) 
(C 22) 
(C 23) 
(C 24) 
(C 25) 
(C 26) 
(C 27) 
(C 28) 
(C 29) 
(C 30) 
(C 31) 
(C 32) 
(C 33) 
(C 34) 
(C 35) 
(C 36) 
(C 37) 
(C 38) 
(C 39) 
(C 40) 
(C 41) 
(C 42) 
(C 43) 
(C 44) 
(C 45) 
(C 46) 
(C 47) 
(C 48) 
(C 49) 
(C 50) 
(C 51) 
(C 52 I 
(C 53) 
(C 54) 
(C 55) 
(C 56) 
(C 57) 
(C 58) 
(C 59) 
(C 60) 
(C 61) 
(C 62) 
(C 63) 
(C 64) 
(C 65) 
(C 66) 
(C 67) 
(C 68) 
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(C 

69) 

(C 

70) 

(C 

71) 

(C 

72) 


YOVERH=( FLOAT (Ml ) -FLOAT (N1 ) ) *S I GMA+GAMMA* (- . 2 ) 

ZQVERH=S I GMA*GAMMA*T AN ( LAMBDA ) *S IN ( ALPHA )*( ABS ( ZSTAR )-ABS ( XS TAR ) ) 

IF ( INK.EQ.7) GO TO 850 
CALL VARLET (C(K)) 

********************* SEE APPENDIX H FOR SUBROUTINE VARLET ********************* 
GO TO 832 (C 73) 

850 CALL DLTAS (C(K)) (C 74) 

********************* see APPENDIX H FOR SUBROUTINE VARLET ********************* 


832 DO 805 Ll=l,28 


(C 75) 


802 

CONTINUE 



(C 77) 

801 

CONTINUE 



(C 78) 


IF ( INK.EQ.7) GO TO 835 



(C 79) 


DELTA(5)=XDELTA( 5 ) * SUML+C0NST1 



(C 80) 


IF ( ABS(DELTA(5) ).LT. 0.00005) GC 

TO 833 


(C 81) 


DEL2 = DEL T A ( 5 ) 



(C 82) 


IF ( INK.NE.l ) GO TO 831 



(C 83) 


DEL1=DEL T A ( 5 ) 



(C 84) 


ZETA1=ZET2=0.9 



(C 85) 


I NK = 2 



(C 86) 


XDEL T A( 5 ) = 0 • 0 



(C 87) 


GO TO 812 



(C 88) 

831 

IF (ABS(DELTA(5) ).LT. 0.00005) GO 

TO 833 


(C 89) 


Z ETA WET A 1- ( ZE T2-Z ET 1 )*DEL2/( DEL2 

-DELI) 


(C 90) 


QEL1 = DEL2 



(C 91) 


ZE T1 = Z ET2 



(C 92) 


I NK=3 



(C 93) 


ZET2=ZETA1 



(C 94) 


XDELT A ( 5 ) = 0 . 0 



(C 95) 


GO TO 812 



(C 96) 

833 

I NK= 7 



(C 97) 


XDELT A ( 5 ) — 0 • 0 



(C 98) 


STORE(K) =Z ET A1 



(C 99) 


GO TO 812 



(C 100) 

835 

DO 807 L3=l t 28 



(C 101) 

807 

DELTA(L3)=XDELTA(L3)*SUML*C0N$T1 



(C 102) 


3ZETA=ZETA1/ < 1 . 0-AC*ZETAl*SI N( ALPHA) ) 


(C 103) 


HCL=( l.O/BZETAJ-l.O 



(C 104) 


WRITE (6,149) C (K) ,BZETA,HCL 



(C 105) 


WRITE (6,150) ( DELTA( I ) , 1=1,25,4) 



(C 106) 


WRITE (6,151) (DELTA(I), 1=2, 26,4) 



(C 107) 


WRITE (6,152) (DELTA! I), 1=3,27,4) 



(C 108) 


WRITE (6,153) (DELTA( I ), 1=4,28,4) 



(C 109) 


IF ( SIGMA. EQ.O. ) GO TO 822 



(C 110) 


IF (C(K) • EQ . 90. ) GO TO 820 



(C 111) 


XS1=TAN(C(K)*RAD )/( 2 . 0*S I GMA*GAMM A*BZ ETA ) 


(C 112) 


CHIM=2.0*C(K)-90. 



(C 113) 


XS2=T AN(CHIM*RAD )/ ( 2 . 0 *S I GMA*GAMMA*BZETA ) 


(C 114) 


WRITE (6,120) X SI , XS2 



(C 115) 


GO TO 822 



(C 116) 

820 

WRITE (6,119) 



(C 117) 

822 

DO 814 L4=l , 28 



(C 118) 

814 

XDELTA( L4 ) =0 • 



(C 119) 

41 

CONTINUE 



(C 120) 


IF ( ITAIL.EQ.O) GO TO 1 



(C 121) 


SUML=0. 031526 



(C 122) 


IF (LI.EQ.l) SUML=0 .025 



(C 123) 


IF ( SIGPAW.EQ.O. . AN C. S I GMAT • EQ.O.) 

GO TO 

450 

(C 124) 


IF (SIGMAW.EQ.O. . AN D . S I GMAT . NE . 0 . ) 

GO TO 

455 

(C 125) 


IF (SIGMAW.NE.G. . AND. S I GMAT . EQ .0 . ) 

GO TO 

460 

(C 126) 


M 7=4 



(C 127) 


APPENDIX C 


N7 = 10 

CONST 1 = 1*0 

IF (ETAl.NE.l.) GO TO 412 
M7=2 

CONST 1 = 2*0 
GO TO 412 
450 M7=N7=1 

XLOAO ( 1 ) = 1 • 0 
SUML=0 .025 
CONST 1 = 40 • 0 
GO TO 412 
455 M7=4 
N7=l 

XLOAO ( 11=1*0 
SUML=0*025 
CONST 1 = 10,0 

IF (ETAl.NE.l.) GO TO 412 
M7=2 

CONST 1=20 • 0 
GO TO 412 
460 M 7 = 1 
N7 = 10 

CONST 1=4 • 0 

IF (ETAl.NE.l.) GO TO 412 
N7 = 5 

CONST 1 = 8.0 
DO 400 L 2=1 » 28 
XOE L T A ( L2 ) =0 . 

400 DELT A ( L2 ) = Q • 

412 DO 420 ITL=1,ITAIL 
TL=ATL*FLOAT ( ITL )+AC 
TH=ATH*FLOAT ( ITL > 

WRITE (6,901) SIGMA, IALPHA,ALAM .GAMMA, ETA1 ,AALP 
WRITE (6,210) 

WRITE (6,211) 

WRITE (6,212) 

WRITE (6,213) 

WRITE (6,214) 

WRITE (6,215) 

WRITE (6,216) 

WRITE (6,217) 

WRITE (6,218) 

DO 91 K= 1,8 
ZETA1=ST0RE(K) 

DO 401 M1=1,M7 

DO 402 N 1= 1 , N7 
XSTAR=(2.*FL0AT (Nl)-U.)/10.0 
YSTAR= ( 1 1 *— 2.*FLQAT (Nl) ) / 10 .0 
ZSTAR=(5.-2.*FLOAT( Ml) )/4.0 
ET A=ETA1+X$TAR*S IGMA 

Z ETA= Z ET A 1 / ( 1 .-A BS ( YST AR ) *SI GMA*GAMMA*ZE T A1*T AN( LAMBDA )* 

1 SIN( ALPHA)) 

XOVERH=TL*COS(ALPHA)+TH*SIN( ALPHA ) -SIGMA *GAMMA*T AN (LAMBDA)* 
1 CUS( ALPHA) *ABS (YSTAR) 

YOVERH=ZSTAR*SIGMAT*GAMMA-YSTAR*SI GMA+GAMMA 

ZGVERH=TH*COS(AL PHA )-TL*SIN< ALPHA ) +S I GMA*GAMMA*T AN ( L AMBD A ) * 
1 S IN ( ALPHA) *ABS (YSTAR ) 

CALL DLTAS (C(K)) 


(C 128) 
(C 129) 
(C 130) 
(C 131) 
(C 132) 
(C 133) 
(C 134) 
(C 135) 
(C 136) 
(C 137) 
(C 138) 
(C 139) 
(C 140) 
(C 141) 
(C 142) 
(C 143) 
(C 144) 
(C 145) 
(C 146) 
(C 147) 
(C 148) 
(C 149) 
(C 150) 
(C 151) 
(C 152) 
(C 153) 
(C 154) 
(C 155) 
(C 156) 
(C 157) 
(C 158) 
(C 159) 
(C 160) 
(C 161) 
(C 162) 
(C 163) 
(C 164) 
(C 165) 
(C 166) 
(C 167) 
(C 168) 
(C 169) 
(C 170) 
( C 171 ) 
(C 172) 
(C 173) 
(C 174) 
(C 175) 
(C 176) 
(C 177) 
(C 178) 
(C 179) 
(C 180) 
(C 181) 
(C 182) 
(C 183) 
(C 184) 
(C 185) 


******** ************* SEE APPENDIX G FOR SUBROUTINE DLTAS ********************** 


DO 405 L 1= 1 , 28 


(C 186) 


405 XDELTA(L1 ) =XDEL TA(L1) +(DELTA(L1 >*XL0AD(N1) ) 


(C 187) 


41 
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402 

CONTINUE 



(C 

188) 

401 

CONTINUE 



(C 

189) 


DO 407 L 3=1 * 28 



(C 

190) 

407 

DELTA(L3)=XDELTA (L3 )*SUML*C0NST1 



(C 

191) 


BZETA=ZETAl/( 1.0 -AC*ZETA1*S IN( ALPHA) ) 



(C 

192) 


TTL=ATL*FLOAT(ITL) 



(C 

193) 


WRITE (6,148) C(K> ,BZETA, S I GMAT , TTL , TH 



(C 

194) 


WRITE (6,150) ( DELTA(I) ,1=1,25,4) 



(C 

195) 


WRITE (6,151) (DELTA! I) ,1=2,26, A) 



(C 

196) 


WRITE (6,152) ( DELTA( I) ,1=3,27,4) 



(C 

197) 


WRITE (6,153) (DELTA(I) ,1=4,28,4) 



<c 

198) 


DO 414 L4=l , 28 



(C 

199) 

414 

XDELTA(L4)=0*0 



(C 

200) 

91 

CONTINUE 



(C 

201) 

420 

CONTINUE 



(C 

202) 


GO TO 1 



IC 

203) 

119 

FORMAT ( 5X*IMPI NGEMENT DISTANCE IS INFINITE*) 



(C 

204) 

120 

FORMAT ( 5X* IMP I NGEMENT OISTANCE (SPANS) USING CHI 

(EFFECTIVE) IS* 

(C 

205) 

1F10.2,5X*USING CHI I MOMENTUM ) IS+F10.2) 



(C 

206) 

148 

FORMAT { // 5H CH I =F5 • 1 , 5X 5HZETA=F7.3, 5X*I NTERFERENCE 

AT 

TAIL+5X 

(C 

207) 

1*SIGMA( TAIL)=*F6*3»5X*TL /H=*F7 . 3 , 5X*TH/H=*F7 . 3 / ) 



(C 

208) 

149 

FORMAT ( // 5H CHI =F5 .1 , 5X 5HZETA=F7.3, 5X*I NTERFERENCE 

AT 

WING* 

(C 

209) 

15X*H(Cl)/H=*F6.3/> 



(C 

210) 

150 

FORMAT (3X5HI W,L )7I F17.4) ) 



(C 

211 ) 

151 

FORMAT (3X5H(U,L )7<F17.4 ) ) 



(C 

212) 

152 

FORMAT (3X5H(W,0)7( F17.4) ) 



(C 

213) 

153 

FORMAT (3X5H(U,0)7(F17.4)) 



(C 

214) 

210 

FORMAT (1X13K1H-H 



(C 

215) 

211 

FORMAT (1X1HI11X1HI31X61HCORRECTICN FACTORS FOR CORRECTING FROM A 

(C 

216) 


1WIND TUNNEL WHICH IS25X1HI) (C 217) 

212 FORMAT ( 1 X 1H 1 11 X 1H 1 1 1 7 ( 1H- ) 1 H I) (C 218) 

213 FORMAT ( 1 X 1H 1 1 1 X 1H 1 16X 1H I 5X6HCL0SED5 X1H I 1 6 X 1H I 2X 1 2HC LOS ED FLOOR2X1 (C 219) 

1HI6X4H0PEN6X1HI 1 6X1HI 5X6HCLOSED4X1HI ) (C 220) 

214 FORMAT ( 1 X 1H I 3X 5 HDELTA3X 1HI 5X6HCL0SED5X1 HI 4X9H0N BOTTOM3X1H I 6X4H0P (C 221) 

1 EN6X1 HI 6X4H0NL Y6 X1H I 5X5HFL00R6XIHI5X6HCL CSE05X1HI3X9H0N 80TT0M3X 1H (C 222) 
21) <C 223) 

215 FORMAT ( 1 X 1H 1 11 X 1H I 16X1H 1 6X4HGNL Y6X 1H 1 16X 1 8H I ( GROUND EFFECT) I6X4H (C 224) 

10NLY6X1HI 16X1HI6X4HCNLY5X1HI ) (C 225) 

216 FORMAT ( 1 X 1 H 1 11 X 1H I 84 ( 1H- ) 1H 1 32< 1H- > 1H I ) (C 226) 

217 FORMAT ( 1 X 1 H I 1 1 X 1H I 36X1 1 HTO FREE A IR37X1 HI 8X16HT0 GROUND EFFECT8X1 (C 227) 

1HI ) (C 228) 

218 FORMAT (1X13K1H-)) <C 229) 

900 FORMAT ( I 1 , F9 * 3 , 4F 1 0 . 3/1 1 , F9 . 3 » 2F 1 0 . 3 ) (C 230) 

901 FORMAT ( 1 H 1 / //40 X* A VER AG E INTERFERENCE OVER A SWEPT WING OF FINITE (C 231) 

1 $PAN*//31X*ZETA COMPUTED TO YIELD ZERO DELTA ( W , L ) IN A CLOSED-ON- (C 232) 
2B0 TTOM-ONL Y TUNNEL*///3GX*SIGMA =*F6 . 3 , 1 5X , A8* L0ADING*15X, (C 233) 

3*L AMBDA =* F7 • 3/ / 30X *GAMM A =* F6. 3 , 17X* ETA =*F7*3,17X* ALPHA = * (C 234) 

4F7*3//41X*ALL DIMENSIONS MEASURED FROM THE AERODYNAMIC CENTER*//) (C 235) 

999 STOP (C 236) 

ENO (C 237) 


42 



APPENDIX D 


FORTRAN PROGRAM FOR COMPUTING INTERFERENCE OVER ROTORS 
IN A VARIABLE MODEL-HEIGHT CLOS ED- ON-BOT TOM-ON LY TUNNEL 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN, VERSION 2.1, TO RUN ON CDC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


INPUT IS AT ADDRESS 1 IN FORMAT 900 (TWO CARDS PER CASE). THE REQUIRED 
INPUT VARIABLES ARE 

LI DISK-LOAD DISTRIBUTION INDICATOR, LI=1 FOR UNIFORM LOADING, 

L I = 2 FOR TRIANGULAR LOADING 

ETA1 DISTANCE FRCM ROTOR CENTER TO RIGHT-HAND WALL DIVIDED 

BY TUNNEL SEMIWIDTH 

GAMMA WIDTH-HEIGHT RATIO OF WIND TUNNEL 

SIGMA RATIO OF ROTOR DIAMETER TO TUNNEL WIDTH 

ALPHA ANGLE OF ATTACK OF ROTOR TIP-PATH PLANE, DEG 

ITAIL NUMBER (0 TO 9) OF TAIL POSITIONS, PROPORTIONATELY INCREASED 

IN TAIL LENGTH AND HEIGHT, AT WHICH INTERFERENCE IS 
REQUIRED. USE SMALLEST LENGTH AND HEIGHT AS INPUT. 

SIGMAT RATIO OF TAIL SPAN TO TUNNEL WIDTH 

TL DISTANCE OF TAIL BEHIND ROTOR CENTER, NONDI MENS IONALI ZED WITH 

RESPECT TO ROTOR RADIUS 

TH HEIGHT OF TAIL ABOVE ROTOR CENTER, NOND I ME N S I ONAL IZ ED WITH 

RESPECT TO ROTOR RADIUS 

ALPHAB ANGLE OF ATTACK OF BODY CARRYING TAIL, DEG 


THIS PROGRAM REQUIRES SUBROUTINE DLTAS (SEE APPENDIX GJ AND SUBROUTINE 
VARLET ( SEE APPENDIX H ). 


SINCE THE TAIL LOCATIONS ARE NCND I ME NS I ON AL I ZE D WITH RESPECT TO ROTOR 
RADIUS, THIS PROGRAM CAN NOT ACCOMMODATE, AND WILL REJECT, TAIL INTERFERENCE 
CALCULATIONS FOR CASES INVOLVING SIGMA = 0, HOWEVER, THE: INTERFERENCE AT THE 
CENTER OF LIFT WILL BE CALCULATED FOR SUCH CASES. TAIL INTERFERENCE FOR SUCH 
CASES CAN BE OBTAINED FROM THE CORRESPONDING PROGRAM FOR WINGS, SINCE, WHEN 
SIGMA IS ZERO, THE REPRESENTATION OF WINGS AND ROTORS IS IDENTICAL. 
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PROGRAM W l NO TUN ( INPUT » OUT PUT f TAP E 5- INPUT ,TAPE6=QUTPUT) (0 l) 

COMMON ZETA, ETA, GAMMA, XOVERH , YOVERH , ZOVE RH , DELTA ( 28 ) (D 2) 

DIMENSION XDELTA(2 8 ),PSI ( 20 > , XLOAD ( 20 > ,RUNIF(20> ,RTRIA( 20) , (D 3) 

1 C ( 8 ) , STORE ( 8 ) (D 4) 

DATA (RUNIF( I), 1=1, 2G)/4*0. 2981, 8*0.6255,8*0.8921/ (D 5) 

DATA (RTRIA(I), 1=1, 20)/4*0.4366,8*0.7296, 8*0.9262/ (D 6) 

DATA (C( I),I=l,8)/20.,30.,40.,50.,60.,70.,80.,90./ (D 7) 

RAD=0. 0174532925199 (D 8) 

PI=3. 14159265358979 (D 9) 

P$I(l)=(PI/4.) (D 10) 

PS I ( 2 )=3. *PSI ( 1 ) (D 11) 

PSI ( 3 )=5 • *PS I ( i ) (D 12) 

PSI(4)=7.*PSI(1> (D 13) 

PSI(5)=PSI(13)=(PI/8.) (D 14) 

P$I(6)=PSI ( 14)=3.*PSI(5) (0 15) 

PSI(7)=PSI(15)=5.*PSI(5) (D 16) 

PSI(8)=PSI(16)=7.*PSi(5) (D 17) 

PSI(9)=PSl(17)=9.*PSI(5) (D 18) 

PSI(10)=PSI(18)=li.*PSI(5) < D 19) 

PS l ( 11 ) = P S I ( 19) = 13.*PSI( 5) (D 20) 

PSI(12)=PSI(20)=15.*PSI(5) (0 21) 

1 READ (5,900) L I ,E TA1 , GAMMA , S I GMA , ALPHA, IT AI L , SI GHAT , TL , TH, ALPHAS (0 22) 

IF (EOF, 5) 999,47 (D 23) 

47 IF (LI.EQ.1) GO TO 804 (D 24) 

I ALPHA=10HTRI ANGULAR (D 25) 

DO 808 M2= 1,20 (D 26) 

808 XLOAD ( M2 ) = RTR I A( M2 ) (0 27) 

GO TO 160 (D 28) 

804 I ALPHA=10H UNIFORM <D 29) 

DO 809 M2=l , 20 (0 30) 

809 XLOAD ( M2 ) =RUN IF ( M2 ) (D 31) 

160 WRITE (6,901) I ALP HA , S I GM A, ALPHA , GAMMA , S I GM AT , ALPHAS , ET Al (D 32) 

WRITE (6,210) (D 33) 

WRITE (6,211) (D 34) 

WRITE (6,212) (D 35) 

WRITE (6,213) (D 36) 

WRITE (6,214) (D 37) 

WRITE (6,215) (D 38) 

WRITE (6,216) (D 39) 

WRITE (6,217) (D 40) 

WRITE (6,218) (D 41) 

SUML=0.0025 (D 42) 

CONS T 1=1 . (D 43) 

ATL=TL (0 44) 

ATH=TH (D 45) 

ALPR= ALPHA (0 46) 

ALPB= ALPHAB (D 47) 

ALPHA=ALPR*RAD (D 48) 

ALPHAB=ALPB*RAD (D 49) 

DO 41 K= 1,8 (D 50) 

DO 803 L 1= 1 , 28 (D 51) 

DELTA(L1)=0. (D 52) 

803 XDEL T A ( L 1 ) =0 • (D 53) 

I NK= 1 (D 54) 

ZETA1=ZET1=1.0 (D 55) 

M7=N7=20 (D 56) 

IF (SIGMA. NE. 0. ) GO TO 811 (0 57) 

N7=M7= 1 (0 58) 

CONST 1=4C0 (D 59) 

GO TO 812 ( D 60 ) 

811 IF (ETA1.NE.1.) GO TO 812 (D 61) 
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CONS T 1=2 . (0 62) 

812 00 801 M1=1,M7 (0 63) 

DO 802 N1=1,N7 (D 64) 

836 IF (ETA1.NE.1.) GO TO 840 <D 65) 

IF (PSKNl ).GT.PI) GO TO 802 (D 66) 

840 ETA=ETAl-( XL0AD(N1 ) *S I GM A *S I N ( PS I ( N1 ) ) ) (0 67) 

ZETA=l./(( l./ZETAl)-(XLOAD(Nl)*SIGMA *S I N ( ALPHA ) *COS ( P S I ( N 1 ) ) *GA MM (D 68) 

1A) ) (D 69) 

XOVERH=S I GMA *GAMMA*COS ( ALPHA ) * C XLOAD ( Ml ) ♦COS ( PS I C HI ) ) -XLGAD ( N 1 ) *C (0 70) 
1 OS ( PS I ( N1 ) ) ) ID 71) 

YOVERH=SIGMA ♦GAMMA* (XL0AD(M1)*SIN(PSI(M1) )-XLOAD(Nl )*S IN( PSI (NI ) ) (D 72) 
1) (0 73) 

ZOVERH=-SIGMA *G AMMA*S IN ( ALPHA )* I XLOA D ( Ml ) *COS ( P S I ( M 1 ) ) -XLOAD! Nl ) * (D 74) 

1 COS (PSI(Nl))) 10 75) 

IF (INK.EQ.7) GO TO 850 (0 76) 

CALL VARLET (C(K)) (0 77) 

********************* SEE APPENDIX H FOR SUBROUTINE VARLET ********************* 
GO TO 832 (D 78) 

850 CALL DLTAS (CCK)) (D 79) 

******** ************* SEE APPENDIX G FOR SUBROUTINE DLTAS ********************** 


832 DO 805 Ll=l,28 

805 XDELTA(L1)=XDELTA(L1)+(DELTA(L1) ) 

802 CONTINUE 
801 CONTINUE 

IF { INK.EQ.7) GO TO 835 

DELTA! 5)=XDELTA( 5) *SUML*CONS T1 

IF ( ABS( DELTA(5) ) .LT. 0.00005) GC TO 833 

DEL2=DELTA(5) 

IF ( INK.NE.I) GO TO 831 
DEL1 = DELT A ( 5 ) 

ZETA1=ZET2=0.7 
I NK=2 

XDELTA! 5 ) =0. 0 
GO TO 812 

831 IF (ABS(DELTA(5) ).LT. 0.00005) GC TO 833 
ZETA1=ZETA1-(ZET2-ZET1 ) *DEL2 / ( DEL 2-DE L 1) 
DEL1=DEL2 
ZET1=ZET 2 
I NK= 3 

ZET2=ZETA 1 
XDELT A ( 5 ) = 0. 0 
GO TO 812 

833 I NK= 7 

XDELT A ( 5 ) =0 • 0 
ST0RE(K)=ZETA1 
HCL=( 1.0/ZETA1)-1.0 
GO TO 812 

835 DO 807 L3=l,28 

807 DELTA ( L3>= XDELT A <L3 )*SUML*C0NST1 
WRITE (6,149) C(K) ,ZETA1,HCL 
WRITE (6,150) (DELTA! I ), 1=1,25,4) 

WRITE (6,151) (DELTA! I ), 1=2,26,4) 

WRITE (6,152) (DELTA! I ), 1 = 3,27,4) 

WRITE (6,153) ( DELTA! I ), 1=4,28,4) 

IF (SIGMA. EQ. 0. ) GO TO 822 
IF ( C (K ) • EQ.90. ) GO TO 820 
XS1=TAN(C( K)*RAO ) / ( 2.0 *S I GMA*GAMM A*ZET A1 ) 
CH I M= 2 * 0*C ( K )— 90 • 

XS2=TAN(CHI M*RAD )/ ( 2 • 0 *S I GMA*GAMM A*Z ETA1 ) 
WRITE (6,120) X SI , XS2 
GO TO 822 


(D 80) 
(D 81) 
(D 82) 
(D 83) 
(D 84) 
(D 85) 
(D 86) 
(D 87) 
(D 88) 
(D 89) 
(D 90) 
(D 91) 
(0 92) 
(0 93) 
(D 94) 
(0 95) 
(D 96) 
(D 97) 
(D 98) 
(D 99) 
(D 100) 
(D 101) 
(D 102) 
(D 103) 
(D 104) 
(D 105) 
(D 106) 
( D 107) 
(D 108) 
(D 109) 
(D 110) 
(D 111) 
(D 112) 
(D 113) 
(D 114) 
(D 115) 
(D 116) 
(D 117) 
(D 118) 
(D 119) 
(D 120) 
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820 

WRITE (6,1191 

(D 121) 

822 

DO 814 L 4= 1 , 2 8 

(D 122) 

814 

XDELTA(L4>=0* 

(D 123) 

41 

CONTINUE 

(D 124) 


IF ( ITAIL.EQ.O) GO TO 1 

(D 125) 


IF ( SIGMA. NE.O. I GO TO 400 

(0 126) 


WRITE (6,901) IALPHA,SIGMA,ALPR, GAMMA, S I GMAT , ALPS , ETA1 

(D 127) 


WRITE (6,904) 

(D 128) 


GO TO 1 

(D 129) 

400 

SUML=0.0125 

(D 130) 


IF (SIGMAT. NE. 0. 1 GO TO 411 

(D 131) 


M7 =1 

(D 132) 


CONS T 1=4.0 

(D 133) 


GO TO 412 

(D 134) 

411 

IF (ETA1.NE.1.) GO TO 413 

(D 135) 


M7=2 

(D 136) 


CONST 1=2 • 0 

(0 137) 


GO TO 412 

(0 138) 

413 

M7 =4 

(D 139) 


CONST 1 = 1 • 0 

(0 140) 

412 

DO 420 ITL=1,ITAIL 

(D 141) 


TL=ATL*FLOAT(ITL ) 

(D 142) 


TH=ATH*FLOAT< ITL ) 

(D 143) 


WRITE (6,901 ) I ALP HA, SIGMA, ALPR, GAMMA, S I GMAT, ALPS, ETA1 

(D 144) 


WRITE (6,210) 

(D 145) 


WRITE (6,211) 

(D 146) 


WRITE (6,212) 

(D 147) 


WRITE (6,213) 

(D 148) 


WRITE (6,214) 

(D 149) 


WRITE (6,215) 

(D 150) 


WRITE (6,216) 

(D 151) 


WRITE (6,217) 

(D 152) 


WRITE (6,218) 

(D 153) 


DO 91 K= 1 , 8 

(D 154) 


ZETA1=ST0RE ( K) 

(D 155) 


DO 401 M1=1,M7 

(D 156) 


DO 402 Nl=l *20 

<D 157) 


ETA=ETAl-( XLOAD( N1 ) *S IGM A*S I N ( PS I ( N1 ) ) ) 

(D 158) 


ZETA=l./( (1./ZETA1)-XL0AD(N1 ) *S I GM A*GAMM A*S IN ( ALPHA ) *COS<PSl(Ni) ) ) 

(D 159) 


X0VERH=SIGMA*GAMMA*( (TL*COS( ALPHAB ))+( TH*S I N ( ALPHAS ) )-(XLGAD(Nl) 

(D 160) 

1 

♦COS ( ALPHA ) *COS (PSI(Nl)) ) ) 

(0 161) 


XM1=FL0AT ( Ml ) 

(D 162) 


YOVERH=SIGMA*GAMMA* (-( (2.*XMl-5. )/4. ) *( S I GMAT/S I GMA ) - ( XLOAD( N1 ) 

(D 163) 

1 

*SIN(PSI(N1) ))) 

(D 164) 


ZOVERH=S IGMA*GAMMA* ( ( TH*COS( ALPHAB ) )-( TL *S I N( ALPHAS ) )+( XLOAO(Nl) 

(D 165) 

1 

♦SIN ( ALPHA) *COS (PSI(Nl))) ) 

(D 166) 


CALL DLTAS (C(K)) 

(D 167) 

********************* SEE APPENDIX G FOR SUBROUTINE DLTAS ********************** 


DO 405 L 1 = 1 , 28 

(0 168) 

405 

XDELTA(L1 )=XDELTA(L1)+DELTA( LI) 

(D 169) 

402 

CONTINUE 

(D 170) 

401 

CONTINUE 

(D 171) 


DO 407 L3= 1,28 

(D 1 72 ) 

407 

DELTA(L3)=XDELTA(L3)*SUML*C0NST1 

(D 173) 


WRITE (6,148) C(K),ZETA1,TL,TH 

(D 174) 


WRITE (6,150) (DELTA( I), 1 = 1,25,4) 

(D 175) 


WRITE (6,151) (DELTAU), 1 = 2,26,4) 

ID 176) 


WRITE (6,152) (DELTA( I ), 1=3,27,4) 

(D 177) 


WRITE (6,153) (DELTA(I), 1 = 4,28,4) 

(D 178) 


DO 414 L4= 1 » 2 8 

(0 179) 

414 

XDELTAIL4 )=0*0 

(D 180) 
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91 

CONTINUE 





420 

CONTINUE 






GO TO 1 





119 

FORMAT ( 

5X* I MP I NGEMENT 

DISTANCE 

IS INFINITE*) 


120 

FORMAT ( 

5X*IMPI NGEMENT 

DISTANCE 

(SPANS) USING CHI 

(EFFECTIVE) IS* ( 


1F10.2»5X* USING CHI (MOMENTUM) IS*F10.2) ( 


148 FORMAT ( // 5X*CH I =* F7 . 3 , 5 X *Z ETA = * F7 • 3 , 5 X* I NTE RF E R ENCE AT TAIL*5X ( 

1*TL/R =*F6#3f5X*TH/R =*F6.3/) ( 

149 FORMAT (//5H CH I =F 5 . 1 # 5X 5 HZ E TA=F 7 . 3 , 5X * I NTERF E RENC E AT ROTOR* ( 

15X*H(CL) /H= *F6* 3 / ) ( 

150 FORMAT ( 3X 5H ( W« L ) 7 ( F 1 7 .4 ) ) ( 

151 FORMAT (3X5H(U,L)7( F17.4) ) ( 

152 FORMAT ( 3 X 5H ( W, D ) 7 ( F 1 7 .4 ) ) ( 

153 FORMAT ( 3 X 5H ( U, D ) 7 ( FI 7 .4 ) ) ( 

210 FORMAT < 1X13K1H-) ) ( 

211 FORMAT ( 1X1HI11X1HI31X61HC0RRECTICN FACTORS FOR CORRECTING FROM A ( 

1 MI NO TUNNEL WHICH IS25X1HI) ( 

212 FORMAT ( 1 X 1 H 1 1 1 X 1 H I 1 1 7 l 1 H- ) 1 H I ) ( 

213 FORMAT { 1 X 1H I 1 1 X 1 H I 1 6X 1H I 5 X6 HC LO S E 05X 1 H 1 1 6X 1H I 2X 1 2HC LOS E D FL00R2X1 ( 

1HI6X4H0PEN6X1HU6X1HI5X6HCLOSED4X1HI ) ( 

214 FORMAT ( 1 X 1H I 3X 5 HD ELT A 3X 1 H I 5 X6HC LC SE D5 XI H 1 4 X 9H0N BOT TOM 3 X 1 H I 6X4H0P ( 

1 EN6X1HI6X4H0NLY6X1H15X5HFL00R6X1HI5X6HCL0SE05X1HI 3X9HON B0TT0M3X1H ( 
21 ) ( 

215 FORMAT < 1 X 1HI 1 1 X 1HI 16X1H I 6X4H0NL Y6X1H 1 16 XI 8H I ( GROUND EFFECT) I6X4H ( 

10NLY6X1HI 16X1HI6X4HCNLY5X1HI ) { 

216 FORMAT ( 1 X 1H I 1 1 X 1H I 84 ( 1H- ) 1H 1 32 ( 1 H- ) 1 H I ) ( 

217 FORMAT ( IX 1H 1 1 1 X 1 H I 36X 1 1 HTO FREE A I R37 XI H I 8X 1 6HT0 GROUND EFFECT8 XI ( 

1 HI ) ( 

218 FORMAT (1X13K1H-)) ( 

900 FORMAT ( 1 1 , F9 .3 , 3F 1 0 . 3/ 1 1 , F9 • 3 ,3F10 . 3 ) ( 

901 FORMAT ( 1 H 1 / / // 4 2X * AV ERAGE INTERFERENCE OVER A ROTOR OF FINITE SPA ( 

1N*//30X*ZETA COMPUTED TO YIELD ZERO OELT A ( W » L ) IN A CLOS ED-ON-BO T T ( 
20M-ONLY TUNNEL* ///55X»A10* DISK L GAD I NG */ /30 X *S I GMA ( ROT OR ) =* ( 

3F6.3? 8X*ALPHA(R0T0R) =*F 7 . 3 f 10 X * GA MMA =*F6 . 3/ / 30 X *S I GMA ( TA IL ) =* ( 

4F6.3,8X*ALPHA(TAIL) =*F 7 . 3 » 1 0 X* E TA =*F6.3//) ( 

904 FORMAT <///40X*TAIL IS OMITTED - EQUATIONS ARE NOT VALID WHEN SIGM ( 
1 A ( ROTOR ) = 0 •*///) ( 

999 STOP ( 

END ( 


D 181) 
D 182) 
D 183) 
D 184) 
D 185) 
D 186) 
D 187) 
0 188) 
D 189) 
0 190) 
D 191) 
D 192) 
D 193) 
D 194) 
D 195) 
D 196) 
D 197) 
D 198) 
D 199) 
D 200) 
D 201) 
D 202) 
D 203) 
D 204) 
D 205) 
D 206) 
D 207) 
D 208) 
D 209) 
D 210) 
D 211) 
D 212) 
D 213) 
D 214) 
D 215) 
D 216) 
D 217) 
D 218) 
D 219) 
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FORTRAN PROGRAM F OR COMPUTING INTERFERENCE OVER WINGS 
IN A VARIABLE WIDTH-HEIGHT-RATIO CLOSED TUNNEL 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN, VERSION 2. It TO RUN ON CDC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


NOTE THAT ALL INPUT AND OUTPUT VARIABLES ARE REFERENCED TO THE GAMMA= 2.0 
TEST SECTION. INPUT IS AT ADDRESS 1 IN FORMAT 900 (TWO CARDS PER CASE). THE 
REQUIRED INPUT VARIABLES ARE 


LI SPAN LQAD-DISTR IBUTICN INDICATOR, LI = 1 FOR UNIFORM LOADING, 

LI =2 FOR ELLIPTIC LCADING 

ZETA1 SEMIHEIGHT OF GAMMA =2 SECTION DIVIDED BY HEIGHT OF AERODYNAMIC 

CENTER ABOVE FLOOR IN GAMMA= 2 SECTION 

ETA1 DISTANCE FROM AERODYNAMIC CENTER TO RIGHT-HAND WALL DIVIDED 

BY TUNNEL SEMIWICTH 

SIGMA RATIO OF WING SPAN TO TUNNEL WIDTH 

LAMBDA WING SWEEP ANGLE, DEG 

ALPHA ANGLE OF ATTACK OF WING, DEG 

ITAIL NUMBER (0 TO 9) OF TAIL POSITIONS, PROPORTIONATELY INCREASED 

IN TAIL LENGTH AND hEIGHT, AT WHICH INTERFERENCE IS 
REQUIRED. USE SMALLEST LENGTH AND HEIGHT AS INPUT. 

SIGMAT RATIO OF TAIL SPAN TO TUNNEL WIDTH 

TL DISTANCE OF TAIL BEHIND AERODYNAMIC CENTER, NOND I MENS IONAL I ZED 

WITH RESPECT TO SEMIFEIGHT OF GAMM A=2 SECTION 

TH HEIGHT OF TAIL ABOVE AERODYNAMIC CENTER, NO ND I M ENS I ONAL I Z E D 

WITH RESPECT TO SEMIHEIGHT OF GAMM A=2 SECTION 


THIS PROGRAM REQUIRES SUBROUTINE DLTAS (SEE APPENDIX G). 


PROGRAM WINDTUNt INPUT, OUTPUT , TAPE 5- INPUT, TAP E6=0UTPUT ) (E 1) 
COMMON ZETA, ETA, GAMMA, XOVERH, YOVERH,ZOVE RH,OEL TA < 28) (E 2) 
DIMENSION XLE (10 ) ,XLQAD( 10),XDELTA(28)»C( 8) (E 3) 
REAL LAMBDA (E 4) 
DATA (CUI, 1 = 1,8 1/20. ,30. ,40. ,50. ,60. ,70. ,80. ,90./ (E 5) 
XLE(1)=XLE(10)=0. 43579 (E 6) 
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XLE(2)=XLE<9)=0. 71422 
XLE(3)=XLE(8)=0. 86603 
XLE(4)=XLE(7)=0. 95394 
XLE(5)=XLE(6)=0. 99499 

X READ ( 5, 900 ) L l , ZETA1 , ET A1 , S IGMA , LAMBDA, AL PH A , I TA IL , S I GMAT , TL , TH 
IF (EOF, 5) 999,47 
47 A ALP= ALPHA 
AL AM=L AMBDA 
ATL=TL 
ATH=TH 
AZETA=ZETA1 
RAD=0. 0174532925199 
ALPHA= ALPHA *R AO 
L AMBDA=LAMBDA*RA C 
00 101 Kl= 1,16 
IF (LI.EQ.1) GO TO 8C4 
IALPHA=8HELLIPTIC 
SUML=0. 0126104 
00 808 M2- 1,10 

808 XL OAO (M2)=XLE(M2) 

GO TO 160 

804 SUML=0.01 

IALPHA=8HUNIF0RM 
DO 809 M2= 1,10 

809 XLOAO ( M2 )= 1 • 0 

160 GAMMA=0.4+0.1*FLOAT(K1) 

AC=0. 5*S I GMA* TAN (LAMBDA) *GAMMA 
IF (LI.NE.l) AC =0 • 848*AC 

WRITE (6, 901 ) I ALPHA, GAMMA, ET A1 , S IGMA , AZ ET A , A AL P , AL AM 
WRITE (6,210) 

WRITE (6,211) 

WRITE (6,212) 

WRITE (6,213) 

WRITE (6,214) 

WRITE (6,215) 

WRITE (6,216) 

WRITE (6,217) 

WRITE (6,218) 

CONST 1 = 1 .0 
DO 41 K= 1 , 8 
DO 803 L 1= 1 , 2 8 
DELTA(L1 )=0. 

803 XDEL TA( L 1 ) =0 • 

IF (SIGMA. NE.O. ) GO TO 811 
M7=N7= 1 
XLOAD(l)sl.O 
SUML= 1 • 

GO TO 812 

811 IF (ETA1.NE.1.) GO TO 813 
M7 = 5 

N7 = 10 
CONST 1 = 2 • 

GO TO 812 

813 M7=N7=10 

812 ZETAl=AZETA/( 1.0 +AC+AZET A*S I N ( ALPHA) ) 

Z ETA 1=2 • 0/ ( 2 • 0+G AMM A* ( 1 • 0— ZE TA1 ) ) 

DO 801 M1=1,M7 

DO 802 N1=1,N7 

XSTAR=(11«— 2«*FLQAT (Ml ) ) /10. 

YSTAR=(2**FL0AT ( N1 ) — 1 1 • ) / 10 • 

ZSTAR=(11.-2.*FL0AT (Nl) ) /10. 


(E 7) 
(E 8) 
(E 9) 
(E 10) 
(E 11) 
(E 12) 
(E 13) 
(E 14) 
(E 15) 
(E 16) 
(E 17) 
(E 18) 
(E 19) 
(E 20) 
(E 21) 
(E 22) 
(E 23) 
(E 24) 
(E 25) 
(E 26) 
(E 27) 
(E 28) 
(E 29) 
(E 30) 
(E 31) 
(E 32) 
(E 33) 
(E 34) 
(E 35) 
(E 36) 
(E 37) 
(E 38) 
(E 39) 
(E 40) 
(E 41) 
(E 42) 
(E 43) 
(E 44) 
(E 45) 
(E 46) 
(E 47) 
(E 48) 
( E 49) 
(E 50) 
(E 51 ) 
(E 52) 
(E 53) 
(E 54) 
(E 55) 
(E 56) 
(E 57) 
(E 58) 
(E 59) 
(E 60) 
(E 61) 
(E 62) 
(E 63) 
(E 64) 
(E 65) 
(E 66) 
(E 67) 
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ETA=ETAi+Y STAR* SIGMA 

(E 68) 


ZETA=ZETAl/( l.-ABS( YSTAR) *S I GMA*GAMMA*ZETA1*TAN( LAMBDA )*S INI ALPHA) 

(E 69) 


1) 

(E 70) 


XOVERH=SIGMA*GAMMA*T AN (LAMBDA )*CCS< ALPHA ) *( ABS ( XSTAR )-ABS( ZSTAR ) ) 

(E 71) 


YGVERH=< FLOAT (Ml )-FLOAT(Nl ) ) *S IGMA*GAMMA*(~. 2) 

(E 72) 


ZQVERH=SI GMA*GAMMA*TAN( L AM8DA) *S IN (ALPHA )*( A8S(ZSTAR )-ABS( XSTAR) ) 

(E 73) 


CALL DLTAS (C(K> ) 

(E 74) 

********************* SEE APPENDIX G FOR SUBROUTINE DLTAS ********************** 


DO 805 L 1= 1 , 28 

(E 75) 

805 

XDELT A(LU=XDELTA( LI ) + (DELTA! LI ) +XLOAD (Nl) ) 

(E 76) 

802 

CONTINUE 

(E 77) 

801 

CONTINUE 

(E 78) 


DO 807 L 3 = 1 ,28 

(E 79) 

807 

DELTA (L3 ) = XDELT A (L3 ) *SUML*C0NST1* ( GAMMA/ 2 • ) 

(E 80) 


WRITE (6,149) C ( K ) 

(E 81) 


WRITE (6,150) (DELTA! I), 1=1, 25,4) 

(E 82) 


WRITE (6,151) ( DELTA( I ), 1=2,26,4) 

(E 83) 


WRITE (6,152) (DELTA(I), 1=3,27,4) 

(E 84) 


WRITE (6,153) (DELTA! I), 1=4,28,4) 

(E 85) 


IF (SIGMA. EQ.O.) GO TO 822 

(E 86) 


IF ( C ( K ) * EQ.90. ) GO TO 820 

(E 87) 


XS1=T AN (COO* RAD )/( 2.0*S I GMA*GAMMA*ZETA1 ) 

(E 88) 


CHiM=2.0*C (K)-90. 

(E 89) 


XS2=TAN(CHIM*RA0 )/( 2. 0*S I GMA*GAMMA*ZETAi ) 

(E 90) 


WRITE (6,120) X SI , XS2 

(E 91) 


GO TO 822 

(E 92) 

820 

WRITE (6,119) 

(E 93) 

822 

DO 814 L4= 1 , 28 

(E 94) 

814 

XDELTAI L4 ) =0 . 

(E 95) 

41 

CONTINUE 

(E 96) 


IF ( I TAIL. EQ.O) GO TO 101 

(E 97) 


SUML=0. 031526 

(E 98) 


IF (LI.EQ.l) SUML=0 .025 

(E 99) 


IF (SIGMAW. EQ.O. .AND. SIGMAT. EQ.O. ) GO TO 450 

(E 100) 


IF (SIGMAW. EQ.O. .AND. SIGMAT. NE.O.) GO TO 455 

(E 101) 


IF (SIGMAW. NE.O. .AND. SIGMAT. EQ.O. ) GO TO 460 

(E 102) 


M7=4 

(E 103) 


N7=10 

(E 104) 


CONST 1=1 .0 

(E 105) 


IF (ETA1.NE.1.) GO TO 412 

(E 106) 


M7=2 

(E 107) 


C0NST1=2 .0 

<E 108) 


GO TO 412 

(E 109) 

450 

M7=N7=1 

(E 110) 


XLOAD ( 1 ) = 1 .0 

(E 111) 


SUML=0.025 

(E 112) 


C0NST1=40 • 0 

(E 113) 


GO TO 412 

(E 114) 

455 

M7=4 

(E 115) 


N7=l 

(E 116) 


XLOADt 1) = 1 .0 

(E 117) 


SUML=0.025 

(E 118) 


CONST 1=10 . 0 

(E 119) 


IF (ETA1.NE.1.) GO TO 412 

(E 120) 


M7=2 

(E 121) 


CONST 1=20 • 0 

(E 122) 


GO TO 412 

(E 123) 

460 

M7=l 

(E 124) 


N7=10 

(E 125) 


C0NST1=4.0 

(E 126) 


IF (ETA1.NE.1.) GO TO 412 

(E 127) 
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[ TAIL 

I ALPHA, GAMMA, ET A 1 , S I GMA , A Z ET A , A AL P , AL AM 


N7 = 5 

CONST 1 = 8 *0 
412 DO 420 ITL=1, 

WRITE (6,901) 

WRITE (6,210) 

WRITE (6,211) 

WRITE (6,212) 

WRITE (6,213) 

WRITE (6,214) 

WRITE (6,215) 

WRITE (6,216) 

WRITE (6,217) 

WRITE (6,218) 

DO 91 K= 1 , 8 

TL = ( A TL* FLOAT (ITL)*GAMMA/2.0)+AC 
TH=ATH*FLOAT( ITL )*GAMMA/2.0 
ZETA1=AZETA/(1.0+AC*AZETA*SIN( ALPHA) ) 

ZETA1=2*0/(2,0+GAMMA*( 1.0-ZETA1) ) 

DO 400 L 2 = 1 , 28 
XDELTAC L2 ) =0 • 

400 DELTA(L2)=0. 

DO 401 M 1= 1 , M7 
DO 402 N 1 = 1 , N7 
XSTAR=(2,*FL0AT (Nll-ll.J/10.0 
YSTAR=(11 m—2 • *FL OAT (N1 ) )/10.0 
ZSTAR=( 5.-2.*FLOAT( Ml) )/4.0 
ETA=ETA1+XSTAR*SIGMA 

Z E TA= Z ET A 1 / ( l.-ABS ( YSTAR)*SI GMA *G AMM A*ZE T A 1 *T AN ( LAMBDA )* 

1 SIN(ALPHA)) 

XOVERH=TL*COS( ALPHA )+TH*SIN( ALPH A ) -S I GMA *G AMMA*T AN ( L AM BD A ) * 
1 cost ALPHA) *ABS( YSTAR) 

YOVERH=ZSTAR*SIGMAT*GAMMA-YSTAR*SI GMA*GAMMA 

ZOVERH=TH*COS( ALPHA )-TL*SIN( AL PH A ) +S I GMA*G AMMA*T AN ( L AM BDA ) * 
1 SIN< ALPHA) *ABS( YSTAR) 

CALL DLTAS (C(K>) 

********************* SEE APPENDIX G FCR SUBROUTINE DLTAS 
DO 405 L 1= 1 , 2 8 

405 XDELTAUl ) = XDELTA(L1) + (DELTA(L1)*XL0AD(N1) ) 

402 CONTINUE 

401 CONTINUE 
DO 407 L3= 1 , 28 

407 DE LTA(L3)=XDELTA(L3 ) *SUML*C0NST1 * ( GAMMA/ 2 .C ) 

TL=ATL* FLOAT! ITL ) 

TH=ATH*FLOAT ( ITL ) 

WRITE (6,148) C (K) ,SIGMAT ,TL,TH 

( DELTA! 1) , 1=1, 25,4) 

( DELTA! I) ,1=2,26,4) 

( DELTA ( I) ,1=3,27,4) 

(DELTA(I) ,1=4,28,4) 


(E 128) 
(E 129) 
(E 130) 

131) 

132) 

133) 

134) 

135) 
(E 136) 
(E 137) 

138 ) 

139) 

140) 

141) 

142) 

143) 

144) 

145) 

146) 

147) 

148) 
( F 149) 

150) 

151) 

152) 

153) 

154) 

155) 

156) 

157) 

158) 

159) 

160) 
161) 
162) 


(E 
( E 
(E 
( E 
(E 


(E 
(E 
(E 
(E 
(E 
(E 
(E 
(E 
( E 
(E 
( E 


( E 
(E 
( E 
(E 
(E 
( E 
( E 
( E 
( E 
( E 
(E 
( E 
( E 


********************** 
(E 163) 


(E 
( E 


WRITE 
WRITE 
WRITE 
WRITE 
DO 414 


(6.150) 

(6.151) 

(6.152 ) 

(6.153) 

L 4= 1,28 


414 XDELTA(L4)=0.0 
91 CONTINUE 
420 CONTINUE 
101 CONTINUE 
GO TO 1 

119 FORMAT (5X*IMPINGEMENT DISTANCE IS INFINITE*) 

120 FORMAT { 5 X*I MP I N GE M ENT DISTANCE (SPANS) USING CHl (EFFECTIVE) IS* 
1F10.2,5X*USING CHI (MOMENTUM) IS*F10.2) 

148 FORMAT (//5H CH I = F5 . 1 , 5X * I NT ERFE R ENC E AT T A I L* 5X*S I GMA ( T A I L ) =* 
1F6.3,5X*TL/H=*F7 . 3 , 5X*TH /H=* F7 • 3/ ) 

149 FORMAT {//5H CH I =F 5 . 1 , 5 X * I NT ERFE R ENC E AT WING*/) 


164) 

165) 
( E 166 ) 

167) 

168) 

169) 

170) 
171 ) 

172) 

173) 

174) 

175) 

176) 

177) 

178) 

179) 
180 ) 
181) 
182) 

183) 

184) 

185) 

186) 
187) 


( E 
( E 
( E 
(E 
( E 
( E 
(E 
( E 
( E 
( E 
( E 
( E 
( E 
(E 
( E 
( E 
( E 
(E 
( E 
( E 
( E 
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150 FORMAT * 3X5H ( W v L ) 7 ( F17.4 ) ) (E 188) 

151 FORMAT (3X5H(U»L)7IF17*4) ) (E 189) 

152 FORMAT (3X5H < M, 0 ) 7< FI 7 .4 ) ) (E 190) 

153 FORMAT ( 3X5H ( U, D ) 7 t F17 .4 ) ) (E 191) 

210 FORMAT UX131UH-)) (E 192) 

211 FORMAT ( 1X1HI 11X 1HI 31X61 HC OR RECTICN FACTORS FOR CORRECTING FROM A (E 193) 

1WIND TUNNEL WHICH IS25X1HI) (E 194) 

212 FORMAT ( 1 X 1H I 11 X 1HI 11 7 ( 1 H- ) 1HI ) (E 195) 

213 FORMAT < 1X1H 1 11 X 1H 1 16X1H I 5X6HCL0SED5X1HI 1 6X 1H I 2X1 2HC LOSED FL00R2X1 (E 196) 

1 H I 6X4HQPEN6X 1HI 16X1HI5X6 HCLO SE04X 1H I ) (E 197) 

214 FORMAT ( 1X1HI3X5H0ELTA3X 1HI 5X6HCLCSED5X1 HI4X9H0N 80T TOM 3X 1 H I 6X4H0P (E 198) 

1EN6X1HI6X4H0NLY6X1HI5X5HFL00R6X1HI 5X6 HCL 0SED5X1HI 3X9 HON B0TT0M3X 1H (E 199) 
2D <E 200) 

215 FORMAT ( 1 X 1HI 11 X 1 H I 16X1H I6X4HC NL Y6X 1H I 16X18H I( GROUND EFFECT) I6X4H (E 201) 

1QNLY6X1HI 1 6X 1HI 6X4HCNL Y5 X1HI ) (E 202) 

216 FORMAT * 1 X 1H 1 11 X 1H I 84 ( 1H- ) 1H I 32( IH- ) 1H I ) (E 203) 

217 FORMAT ( 1 X 1H 1 11 X 1H I 36X1 1 HTO FREE A I R37X1 H I 8X 1 6HT0 GROUND EFFECT8X1 (E 204) 

1HI ) < E 205) 

218 FORMAT (1X13K1H-)) (E 206) 

900 FORMAT C 1 1 ,F9.3, 4F10.3/I 1 ,F9 . 3 1 2F 10. 3 ) (E 207) 

901 FORMAT < 1H 1///26X* AVERAGE INTERFERENCE FOR SWEPT WING IN A VARIA8L (E 208) 

IE WIDTH-HEIGHT-RATIO CLOSED TUNNEL*// 58XA 8* LOADING*// (E 209) 

2 36X*GAMMA =*F6. 3 f 10 X*E TA =*F7. 3, 10X*SI GMA =*F7.3// (E 210) 

336X*ZETA =*F6. 3 . 10X*ALPHA =*F7 . 3 , 10X*L AM8DA =*F7.3// (E 211) 

435X*ALL VALUES ARE REFERENCED TO THE HIGH-SPEED ( GAMMA=2 .0 ) SECTIO (E 212) 
5N*/39X* AL L DIMENSIONS ARE MEASURED FROM THE AERODYNAMIC CENTER*//) <E 213) 

999 STOP (E 214) 

end (E 215) 
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FORTRAN PROGRAM FOR COMPUTING INTERFERENCE OVER ROTORS 
IN A VARIABLE WI D TH-HE I GHT-RAT IC CLOSED TUNNEL 


THIS PROGRAM WAS WRITTEN IN CDC FORTRAN, VERSION 2.1, TO RUN ON CDC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


NOTE THAT ALL INPUT AND OUTPUT VARIABLES ARE REFERENCED TO THE GAMMA= 2*0 
TEST SECTION. INPUT IS AT ADDRESS 1 IN FORMAT 900 I TWO CARDS PER CASE). THE 
REQUIRED INPUT VARIABLES ARE 


LI OISK-LOAD DISTRIBUTION INDICATOR, LI=1 FOR UNIFORM LOADING, 

L I =2 FOR TRIANGULAR LOADING 

ZETA1 SEMIHEIGHT OF GAMMA=2 SECTION DIVIDED BY HEIGHT OF ROTOR 

CENTER ABOVE FLOOR IN G AMMA-=2 SECTION 

ETA1 DISTANCE FRGM ROTOR CENTER TO RIGHT-HAND WALL DIVIDED 

BY TUNNEL SEMIWIDTH 

SIGMA RATIO OF ROTOR DIAMETER TO TUNNEL WIDTH 

ALPHA ANGLE OF ATTACK OF ROTOR TIP-PATH PLANE, DEG 

ITAIL NUMBER (0 TO 9) OF TAIL POSITIONS, PROPORTIONATELY INCREASED 

IN TAIL LENGTH AND HEIGHT, AT WHICH INTERFERENCE IS 
REQUIRED. USE SMALLEST LENGTH AND HEIGHT AS INPUT. 

SIGMAT RATIO OF TAIL SPAN TO TUNNEL WIDTH 

TL DISTANCE OF TAIL BEHIND ROTOR CENTER, NONDI MENS IONAL I ZE D WITH 

RESPECT TO ROTOR RADIUS 

TH HEIGHT OF TAIL ABOVE ROTOR CENTER, NONO l MENS l ONAL I Z ED WITH 

RESPECT TO ROTOR RADIUS 

ALPHAS ANGLE OF ATTACK OF BODY CARRYING TAIL, DEG 


THIS PROGRAM REQUIRES SUBROUTINE CL T AS (SEE APPENDIX Gl. 


SINCE THE TAIL LOCATICNS ARE NCND I MENS I QNAL I ZED WITH RESPECT TO ROTOR 
RADIUS, THIS PROGRAM CAN NOT ACCOMMODATE, AND WILL REJECT, TAIL INTERFERENCE 
CALCULATIONS FOR CASES INVOLVING SIGMA = 0, HOWEVER, THE INTERFERENCE AT THE 
CENTER OF LIFT WILL BE CALCULATED FOR SUCH CASES. TAIL INTERFERENCE FOR SUCH 
CASES CAN BE OBTAINED FROM THE CORRESPONDING PROGRAM FOR WINGS, SINCE, WHEN 
SIGMA IS ZERO, THE REPRESENTATION OF WINGS AND ROTORS IS IDENTICAL. 


APPENDIX F 


PROGRAM WINDTUNIINPLT, OUTPUT ,TAPE5=INPUT *T APE6=0UTPUT ) CF 1) 

COMMON ZETA»ETA, GAMMA , XOVERH.YOVERH, ZOVERH, DELTA! 28) (F 2) 

DIMENSION XDELTA(28),PSI(20),XLOAD(20),RUNIF(20),RTRIA(20), (F 3) 

1 C(8) (F 4) 

DATA (RUNIF(l), 1=1, 20)/4*0. 2981, 8*0. 6255, 8*0. 8921/ IF 5) 

DATA (RTRIAI 11,1=1,201/4*0.4386,8*0.7296,8*0.9262/ (F 6) 

DATA (C(I), 1=1, 81/20. ,30. ,40. ,50. ,60. ,70. ,80. ,90./ (F 7) 

RAD=0. 0174532925199 CF 8> 

PI=3. 14159265358979 (F 9) 

PSI 1 1 ) = ( P I /4. ) IF 10) 

PS 1 1 2 )=3 . *PSI II) IF 11) 

PS 1 1 3 ) = 5 . *PS I 1 1 ) IF 12) 

PS 1 1 4 )=7 . *PS II 1 ) IF 13) 

PSII 5 ) = PSI 1 131=1 PI/8. ) IF 14) 

PSII6)=PSII14)=3.*PSI(5) IF 15) 

PSI17)=PSII15)=5.*PSI15) IF 16) 

PSI(8)=PSII16)=7.*PSl(5) IF 17) 

PSI19)=PSII17)=9.*PSI(5) IF 18) 

PSI!10)=PSI! 181=11. *PSI<5) IF 19) 

PSI111)=PSII19)=13.*PSII5) IF 20) 

PSIU2) = PSI (201 = 15. *PSI15) IF 21) 

1 READ (5,900) L I ,ZE TA1 , ETA1 , SI GMA , ALPH A, I T AI L ,St GMAT ,TL, TH, ALPHAB IF 22) 

IF (EOF, 5) 999,47 IF 23) 

47 AZETA=ZETA1 (F 24) 

ATL=TL IF 25) 

ATH=TH IF 26) 

ALPR= ALPHA IF 27) 

ALPB= ALPHAB (F 281 

ALPHA=ALPR *RAD IF 29) 

ALPHAB=ALPB*RAD (F 30) 

IF (LI.EQ.l) GO TO 804 IF 31) 

IALPHA=10HTRI ANGULAR IF 32) 

DO 808 M2= 1 » 20 IF 33) 

808 XLOAD I M2 ) = RTR I A I M2 I IF 34) 

GO TO 160 IF 35) 

804 I ALPHA=10 H UNIFORM IF 36) 

DO 809 M2= 1,20 IF 37) 

809 XL0A0IM2) = RUN IF I M2 ) IF 38) 

160 DO 421 K 1 = 1 » 16 IF 39) 

GAMMA=0 .4+0. 1*FL0AT(K1) IF 40) 

WRITE (6,901) IALP HA, GAMMA, SIGMA, ALPR,AZETA,SIGMAT,ALP8,ETA1 IF 41) 

WRITE 16,210) IF 42) 

WRITE (6,211) IF 43) 

WRITE (6,212) IF 44) 

WRITE (6,213) <F 45) 

WRITE (6,214) (F 46) 

WRITE (6,215) IF 47) 

WRITE 16,216) IF 48) 

WRITE 16,217) IF 49) 

WRITE (6,218) IF 50) 

SUML=0 • 0025 IF 51) 

CONST 1=1 . IF 52) 

DO 41 K= 1 » 8 IF 53) 

DO 803 L 1= 1 , 28 IF 54) 

DE LTA I L 1 ) = 0. IF 55) 

803 XDELTAt L 1 ) =0. <F 56) 

M7=N7=20 IF 57) 

IF ( SIGMA. NE.O. ) GO TO 811 IF 58) 

N7=M7=1 IF 59) 
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CONST 1=400 


(F 

60) 


GO TO 812 


(F 

61 ) 

811 

IF (ETAl.NE.l.) GO TO 812 


(F 

62) 


CONST 1=2 • 


(F 

63) 

812 

DO 801 M 1 = 1 » M7 


(F 

64) 


DO 802 N1=1,N7 


(F 

65) 

836 

IF (ETAl.NE.l.) GO TO 840 


(F 

66) 


IF (PSI(Nl).GT.P I) GO TO 802 


(F 

67) 

840 

ETA=ETA1-(XL0AD( N1 )*SIGMA *S I N ( PS I ( N1 ) ) ) 


(F 

68) 


ZETAl=2.0/(2.0+GAMMA*( 1.0-AZETA) ) 


(F 

69) 


ZETA=l./< ( 1./ZETA1 >-(XLOAD(Nl )*SIGMA *S I N ( ALPHA ) *COS ( PS I ( N1 J ) *GA MM 

(F 

70) 

1A) I 


(F 

71) 


XOVERH=S I GMA *GAMMA*COS( ALPHA ) * ( XLOAD< Ml ) *COS ( PS I ( HI ) l-XLOAD ( N1 ) *C 

(F 

72) 

1 OS ( PSI ( N1 ) ) ) 


(F 

73) 


YO VERH=S I GMA *GAMMA* ( XLO A0( M 1 ) *S I N « PS l ( M 1 ) )-XLOAD(Nl ) *SI N ( PS I ( N1 )) 

(F 

74) 

1) 


(F 

75) 


ZOVERH=-SIGMA *G AMMA*S IN ( ALPHA 1* ( XLOACl Ml ) *COS (PS I ( M 1 1 J-XL0AD(N1 )* 

(F 

76) 

1CGS ( PS I ( N1 ) ) ) 


(F 

77) 


CALL DLTAS (C(K)) 


(F 

78) 

********************* SEE APPENOIX G FOR SUBROUTINE DLTAS 

+*****+*************** 


DO 805 Ll=l , 28 


(F 

79) 

805 

XDELTAILl ) = XDEL T A(L1)+(DELTA( LI) ) 


(F 

80) 

802 

CONTINUE 


(F 

81) 

801 

CONTINUE 


(F 

82) 

835 

DO 807 L3=l,28 


(F 

83) 

807 

DELTA(L3)=XDELTA!L3)*SUML*C0NSTl*(GAMMA/2.) 


(F 

84) 


WRITE (6*149) C( K) 


(F 

85) 


WRITE (6,150) ( DELTA( I) , 1=1,25,4) 


(F 

86) 


WRITE (6,151) (DELTA( II, 1 = 2,26,4) 


(F 

87) 


WRITE (6,152) (DELTA! I ), 1=3,27,4) 


(F 

88) 


WRITE (6,153) (DELTA! I ), 1 = 4,28,4) 


(F 

89) 


IF (SIGMA. EQ.O.) GO TO 822 


(F 

90) 


IF (C(K) .EQ.90. ) GO TO 820 


(F 

91 ) 


XS1=TAN(C(K)*RAD ) / ( 2 . 0*S I GMA*G AMM A*ZET A1 ) 


(F 

92) 


CHIM=2.0*C ( K )— 90 . 


(F 

93) 


XS2=T AN (CHI M*RAD )/ ( 2.0*S IGMA*GAMKA*ZETA1 ) 


(F 

94) 


WRITE (6,120) X SI , XS2 


(F 

95) 


GO TO 822 


(F 

96) 

8 20 

WRITE (6,119) 


(F 

97) 

8 22 

DO 814 L 4= 1,28 


(F 

98) 

814 

XDE LT A ( L4 ) =0 • 


(F 

99) 

41 

CONTINUE 


( F 

100) 


IF ( ITAIL.EQ.O) GO TO 1 


( F 

101) 


IF ( SIGMA. NE.O. ) GO TO 400 


(F 

102) 


WRITE (6,901) I ALPHA, GAMMA, SIGMA, AL PR , AZ ET A , S I GM AT , 

ALPB , ET A 1 

(F 

103) 


WRITE (6,904) 


(F 

104) 


GO TO 1 


(F 

105) 

400 

SUML=0. 0 125 


(F 

106) 


IF (SIGMAT.NE.O. ) GO TO 411 


(F 

107) 


M7 = l 


(F 

108) 


CONST 1 = 4 .0 


(F 

109) 


GO TO 412 


(F 

110) 

411 

IF (ETAl.NE.l.) GO TO 413 


(F 

111) 


M7=2 


(F 

112) 


CONST 1=2.0 


(F 

113) 


GO TO 412 


(F 

114) 

413 

M7=4 


(F 

115) 


CONST 1=1.0 


(F 

116) 

412 

DO 420 I TL=1 , I T A I L 


(F 

1171 


TL=ATL*FLOAT( ITL) 


(F 

118) 


TH=ATH*FLOAT( ITL ) 


( F 

119) 
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WR 

ITE 

(6,901 ) I ALPHA, GAMMA, SI GMA, ALPR, AZ ETA , S IGMAT , ALPS, ETA1 

(F 

120) 

WR 

ITE 

(6,210) 

(F 

121 ) 

WR 

ITE 

(6,211) 

(F 

122) 

WR 

ITE 

(6,212) 

(F 

123) 

WR 

ITE 

(6,213) 

( F 

124) 

WR 

ITE 

(6,214) 

(F 

125) 

WR 

ITE 

(6,215) 

(F 

126) 

WR 

ITE 

(6,216) 

(F 

127) 

WR 

ITE 

(6,217) 

(F 

128) 

WR 

ITE 

(6,218) 

(F 

129) 

DO 

91 

K= 1 , 8 

(F 

130) 

ZETA 1= 

2.0/(2.0+GAMMA*(1.0-AZETA) ) 

(F 

131) 

DO 

401 

M 1= 1 »M7 

(F 

132) 

DO 

40 2 

Nl = l , 20 

(F 

133) 


399 


ET A=ETA1- ( XLOAD ( N1 )*SIGMA*SIN(PSI(N1) ) ) (F 134) 

ZETA=l*/{ { 1./ZETA1 ) -XLO AD ( N1 ) *S I GMA*GAMMA*S IN ( ALPHA ) *COS { PS I ( N 1 ) ) » ( F 135) 

XOVERH=SIGMA*GAMMA*(( TL*COS ( ALPHAB >)+ t TH *S I N ( ALPHAS ) )-( XLOAD (Nil (F 136) 

1 *COS (ALPHA)* CO SCPSI(Nl)) ) ) (F 137) 

XM1=FL0AT(M1) (F 138) 

YGVERH=S I GMA*GAM MA* (-( ( 2 . *XM 1- 5, ) /4. ) *< S IGHA T/S l GMA ) — < XLOAD< N1 ) (F 139) 

1 * SI N ( P S I (Nl) ))) (F 140) 

ZOVERH=SIGMA*GAMMA*( (TH*COS( ALPHAB) )-( TL*SIN (ALPHAB) )+(XLOAD(Nl) (F 141) 
1 *SIN(ALPHA) *COS ( PS I ( N 1 ) ) ) ) ( F 142) 

CALL DLTAS (C(K)) (F 143) 

* ** * **** ******* ****** see APPENOIX G FOR SUBROUTINE DLTAS ********************** 
DO 405 L 1=1 » 28 

405 XDELTAIL1 )=XDELTA(L1 )+DELTA( LI ) 

402 CONTINUE 
401 CONTINUE 

DO 407 L3= 1 » 28 

407 DELTA* L3 )=XDELTA ( L3 ) *SUML*CONST 1* (GAMMA/ 2 . 0 ) 

WRITE (6,1 48 ) C ( K ) , TL,TH 

(DELTA! I), 1=1,25,4) 

(DELTA(I), 1=2,26,4) 

(DELTAt I), 1=3,27,4) 

( D ELTA( I), 1=4,28,4) 


WRITE 
WRITE 
WRITE 
WRITE 
DO 414 


(6.150) 

(6.151 ) 

(6.152) 

(6.153) 
L4= 1 ,28 


414 

91 

420 

421 


119 

120 


150 

151 

152 

153 
210 
211 


XDELT A( L4 ) = 
CONTINUE 
CONTINUE 
CONTINUE 
GO TO 1 
FORMAT ( 
FORMAT ( 


= 0.0 


IS INFINITE*) 

( SPANS) USING CHI 


(EFFECTIVE) IS* 


5X*I MPl NGEMENT DISTANCE 
5X* I MP I NGEMENT DISTANCE 
1F10.2,5X*USING CHI (MOMENTUM) IS+F10.2) 

148 FORMAT t//5X*CHI =* F7 .3 , 5X*I NTERFERENCE AT TAIL*5X*TL/R =*F9.3, 

1 5X+TH/R =*F8 • 3/ ) 

149 FORMAT {//5H CHI =F7 . 3 , 5X* I NT ERFE RENCE AT ROTOR*/) 

FORMAT (3X5H(W,L)7(F17.4) ) 

FORMAT (3X5H(U,L)7(F17#4) ) 

FORMAT (3X5H( W,D)7(F17.4) ) 

FORMAT (3X5H(U,D)7(F17.4) ) 

FORMAT ( 1X131(1 H— ) I 

FORMAT ( 1X1HI11X1HI31X61 HCORREC T I CN FACTORS FOR CORRECTING FROM A 
1WIND TUNNEL WHICH IS25X1HI) 

212 FORMAT ( 1 X 1H 1 11X 1H 1 11 7 ( 1 H- )1HI ) 

213 FORMAT ( 1 X 1H I 1 1 X 1H I 16 X 1H I 5X6HCL0S ED5X1 H I 16X 1H I 2X 1 2HCL0S ED FL00R2X1 
1HI6X4H0PEN6X1HI 16X1HI 5X6HCL0SED4X 1H I ) 

214 FORMAT ( 1 X 1H I 3X 5 HD ELT A3X 1H I 5 X6HCLCS ED5X1 HI 4X9H0N B0TT0M3 XI H I 6X4H0P 
1EN6X1HI6 X4H0NLY6 X1H I5X 5HFLOOR6X 1 H I 5X6HCL0SED5X1HI 3X9 HON B0TT0M3X1H 

21 ) 


( F 

144) 

(F 

145) 

(F 

146) 

( F 

147) 

(F 

148) 

(F 

149) 

( F 

150) 

(F 

151) 

(F 

152) 

( F 

153) 

(F 

154) 

(F 

155) 

(F 

156) 

(F 

157) 

(F 

158) 

(F 

159) 

(F 

160) 

(F 

161) 

(F 

162) 

(F 

163) 

(F 

164) 

(F 

165) 

(F 

166) 

(F 

167) 

( F 

168) 

(F 

169) 

(F 

170) 

(F 

171) 

(F 

172) 

(F 

173) 

(F 

174) 

(F 

175) 

(F 

176) 

(F 

177) 

(F 

178) 

(F 

179) 
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215 FORMAT (1X1HI11X1HI 16X1HI6X4H0NLY6X1HI 16X1 8H l ( GROUND EFFECT) I6X4H (F 180) 

10NLY6X1HI 16X1HI 6 X4HCNLY5 X1HI ) (F 181) 

216 FORMAT ( 1 X1H i 11 X 1H I 84 ( 1H- ) 1HI 32C 1H- ) 1 H I ) (F 182) 

217 FORMAT < 1 X 1H 1 1 1 X 1H I 36X1 1 HTO FREE A I R 37 XI H I 8X 1 6HT0 GROUND EFFECT8X1 (F 183) 

1HI) (F 184) 

218 FORMAT (1X131I1H-)) (F 185) 

900 FORMAT ( 1 1 , F9 .3 f 3F 1 0 . 3/ 1 1 , F9 . 3 , 3F 10 . 3 ) (F 186) 

901 FORMAT ( 1H1 /42X* AV ERAGE INTERFERENCE OVER A ROTOR OF FINITE SPA <F 187) 

IN*// 43X* I N A VARIABLE W I DTH- HE IGHT-R AT 1 0 CLOSED T UNNEL*/ // 42 X , A1 0 (F 188) 
2* DISK LOAD I NG*1 OX* GAMMA =*F6 . 3/ /30X*S I G MA ( ROTOR ) =* (F 189) 

3F6.3, 8X*ALPHA(R0T0R ) =*F 7 . 3, 10 X*ZETA = *F 6. 3 // 30X*S I GMA ( T A I L > =* (F 190) 

4F6.3t8X*ALPHA(TAIL) =*F 7 . 3 , 10X* E TA =*F6 . 3 // 3 5X *AL L VALUES ARE RE (F 191) 
5FERENCED TO THE HIGH-SPEED <GAMMA=2.0) SECTION*//) (F 192) 

904 FORMAT (///40X*TAIL IS OMITTED - EQUATIONS ARE NOT VALID WHEN SIGM (F 193) 
1 A( ROTOR ) = 0.*///) (F 194) 

999 STOP (F 195) 

END (F 196) 
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SUBROUTINE DLTAS 


THIS SUBROUTINE WAS WRITTEN IN CDC FORTRAN, VERSION 2.1* TO RUN ON CDC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS SUBROUTINE 
HAS BEEN FOUND TO BE SATISFACTORY CN THE AFOREMENTIONED COMPUTERS WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


THIS SUBROUTINE IS REQUIRED BY THE PROGRAMS OF ALL PRECEEDING APPENDICES. 


SUBROUTINE OLTAS ( ANGL ) 

COMMON ZETA,ETA, GAMMA, XOVERH , YOVERH , ZOVERH, DELTA (28) 
DIMENSION V ( 3 * 9 ) , ADEL ( 28 ) 

SC=SIN< ANGL*Q. 01 74532925199) 

CC=COS( AN GL*0.01745 32 925199) 

Z6=ZETA*Z0VERH+1 . 

Z8=-Z6 

Z7=Z8-1. 

DO 8 J 1=1 , 28 
8 DELTA ( J1 ) =0 . 

DO 10 M= 1 , 7 
DO 10 N= 1,7 

IF (N.EQ.4.AND.M.EQ.4) GO TO 10 
DO II J 1 = 1 * 3 
DO 11 J2=l , 9 

11 V( J1 , J2 ) =0 • 

DO 12 Jl=l , 28 

12 ADEL ( J 1 ) =0 • 

AM=M— 4 
AN=N-4 

X=ZETA*XOVERH 

Y=ZETA*( YOVERH-2 . *AM*GAMMA+GAMMA* ( 1 . - ET A ) * ( 1 - 1 . )**M) ) 
Z=ZETA*( Z0VERH-4.*AN) 

A=SQRTCX*X+Y*Y+Z*Z) 

B = A+Z *CC-X*SC 

V( 1,1 ) = ( ( X *X + Y* Y )/( B*A*A*A) )-( ( Z +A*CC ) / ( B*A ) )**2 
VC 2,1 )=-( X*Z )/C B*A*A*A)-( Z+A*CC) *( X-A*SC )/(B*B*A*A) 

VC 3,1 ) = ( { Y*Y+Z*Z )/( B*A*A*A))-{ (X-A*SC )/C B*A > )**2 
Z=-Z-2. 

A=SQRT(X*X+Y*Y+Z*Z) 

B=A+Z*CC-X*SC 

VC 1,3)= ( { X *X + Y* Y 1/ C B*A*A*A ) )- ( C Z+A*CC ) /< B*A) )**2 
V( 2,3) =-< X*Z )/C B *A* A*A ) - ( Z + A*CC ) * ( X-A*SC )/(B*B*A*A) 

VC 3 , 3 ) = ( <Y*Y + Z*Z )/C B*A*A*A) ) - ( { X-A*SC ) /< B*A ) )**2 

IF (ANGL. EQ. 90.0) GO TO 13 

X=X-($C/CC) 

Z»-Z-l. 

A=SQRT(X*X+Y*Y+Z*Z) 

B=A+Z*CC-X*SC 

V( 1,2 ) = ( C X*X + Y*Y )/( B*A*A*A) )-( (Z + A*CC)/< B*A) )**2 


CG 1) 
(G 2) 
(G 3) 
CG 4) 
CG 5) 
(G 6) 
(G 7) 
CG 8) 
CG 9) 
(G 10) 
CG II ) 
CG 12) 
(G 13) 
(G 14) 
CG 15) 
CG 16) 
(G 17) 
(G 18) 
CG 19) 
(G 20) 
(G 21) 
(G 22) 
(G 23) 
(G 24) 
CG 25) 
(G 26) 
(G 27) 
(G 28) 
CG 29) 
(G 30) 
(G 31 ) 
CG 32) 
CG 33) 
CG 34) 
CG 35) 
(G 36) 
CG 37) 
(G 38) 
(G 39) 
CG 40) 
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V(2f2)=-(X*Z)/(B*A*A*A)-(Z + A*CC)*(X-A*SC)/!B*B*A*A) 

V ( 3 » 2 ) = ( ( Y*Y+Z*Z )/( B*A*A*A) )-( !X-A*SC )/! B*A) )**2 
B= A— X 

V ( 1 1 5 ) = ( <X*X+Y*Y )/( B * A*A *A))—!Z/(B*A) ) ** 2 
V(3*5)=X/( A*A*A) 

Z=-Z 

B=A+Z*CC-X*SC 
V( 2*5)=Z/(A*A*A) 

V! 1,4)=( ! X*X+Y*Y )/( B*A*A*A) )-! (Z+A*CC )/( B*A ) ) **2 

V ( 2 ? 4 ) =— (X*Z)/(B *A* A*A)-(Z+A*CC) * ( X— A* SC )/(B*B*A*A) 

V! 3t4)=( ( Y*Y+Z*Z )/( B*A*A*A) )-( (X-A*SC)/! B*A) )**2 

13 ADEL(l)=V(l f l)-V(l f 2)-V< 1,3) +V! 1,4) 

ADEL!2)=V(2,1)-V (2 * 2 ) +V< 2, 3>-VC 2 ,4 ) 

ADEL!3l=V! 2,1)-V(2,2)-V( 2,3) +V( 2 , 4 ) +2 . *V ! 2 , 5 ) 

ADEL(4)=V(3, 1)-V!3,2)+V( 3,3)-V( 3 , 4 ) +2 . *V ( 3 , 5) 

ADEL (5)=( !-l.)**(M+N) >*ADEL< 1) 

ADEL(6) = ( (-1. )**CM+N) >*ADEL( 2) 

ADEL(7 ) = ( ! — 1 • ) ** ! M + N ) )*AOEL( 3) 

ADEL ( 8 ) = ( ( — 1 • !♦* ! M+N ) ) *AOEL ! 4 ) 

ADEL ! 9 ) = C (-1. )**M)*!V(1,1)-V( i,2)+V(l,3)-V( 1,4)+2.*VC 1» 511 
ADEL ( 10)= ( !-l.)**M)*( V(2,l)-V(2,2>-V(2,3)+V(2t4)+2.*V!2,5I ) 
ADEL ( 11 )=( (-1. ) **M)*I V(2, 1 )-V(2,2)+V(2,3)-V!2,4) ) 

ADEL ( 12)= ( !-l.)**M)*(V(3,l)-V(3,2)-V!3,3>+V|3,4) I 
DO 14 J 1= 1 , 12 

14 DELTA! J1)=DELTA(J1)+ADEL(J1) 

10 CONTINUE 

DO 15 J 1 = 1 , 8 

15 DELTA! J 1+20 ) =DE L TA ( J1 ) 

X= ZET A*XOVERH 
Y=ZETA*YOVERH 

Z = Z7 

A=SQRT!X*X+Y*Y+Z*Z) 

B= A+Z *CC — X*SC 

VC1»7) = (( X*X+Y*Y )/( B*A*A*A) )-( (Z+A*CC )/< B*A) )**2 
V(2,7)=-(X*Z>/( B*A*A*A)-(Z+A*CC)*(X-A*SC )/(B*B*A*A) 

VC 3,7) = ! ( Y*Y+Z*Z )/( B*A*A*A ) ) - ! ( X-A*SC ) / ( B*A ) )**2 

IF ( ANGL.EQ.90.0) GO TO 16 

X=X-(SC/CC) 

Z = Z6 

A= SQR T ( X*X+Y*Y+Z *Z ) 

B=A+Z*CC-X*SC 

V! 1 ,6)=( ( X*X+Y*Y )/{ B*A*A*A ) )-( ( Z+A+CC )/( B*A ) ) **2 
V( 2,6)=-! X*Z)/( B*A*A*A)— ( Z + A*C C ) *( X— A* SC )/(B*B*A*A) 

V( 3,6 ) = ( ( Y*Y+Z*Z )/ ( B*A*A* A ) )-( ( X-A+SC )/( B*A ) )**2 
B= A- X 

V( 1,9 )=( ( X* X+ Y* Y )/( B*A*A*A) ) -!Z/ t B*A ) ) ** 2 
V!2,9)=Z/( A*A*A) 

V!3,9)=X/! A*A*A) 

Z = Z8 

B=A+Z*CC-X*SC 

V! 1,8)=! ! X * X + Y*Y )/ ( B*A*A*A ) )-( ( Z+A*CC )/( B*A) )**2 
V( 2t 8 )=-(X*Z )/( B*A*A*A)-(Z+A*CC)*! X-A*SC )/(B*B*A*A) 

VC 3»8 ) = ! ( Y*Y+Z*Z )/( B*A*A*A) )-( (X-A*SC)/( B*A) 1**2 

16 DELTA! 13 )=-V ! 1 , 6 )-V ! l f 7) +Vt It 8) 

DELTA ( 14)=-V!2,6)+V(2, 7)-V(2 ,8) 

DELTA! 15 )=-V(2,6 )-V(2,7) +V ( 2 , 8 ) +2 . *V ! 2 ,9 ) 

DELTA (16 )=-V!3,6)+V(3,7)-V(3,8)+2.*V(3,9) 

DELTA! 17) = -V( 1,6) +V( 1,7) -V! 1,8 ) + 2.*V< 1,9) 

DELTA! 18)=DELTA< 15) 

DELTA! 19 )=DELTA( 14) 

DELTA! 20)=-V(3,6)-V(3,7) +V(3 f 8) 

DO 17 J 1= 1 » 4 

17 DELTA! J1)=0ELTA( J1 ) +CELT A < J 1 + 1 2 ) 

DO 18 Jl= 5, 12 

18 DELTA! J1 )=DELTA( J1 ) + DELT A! J1 +8 ) 

AMT=— 2« *GAMMA*Z E TA*ZETA/3* 14159265 35 8979 
DO 19 J 1 = 1 1 2 8 

19 DELTA! J1)=AMT*DELTA!J1) 

RETURN 

END 


(G 41) 
(G 42) 
( G 43) 
(G 44) 
(G 46) 
(G 47) 
(G 48) 
(G 45) 
(G 49) 
(G 50) 
(G 51) 
(G 52) 
! G 53) 
(G 54) 
(G 55) 
(G 56) 
! G 57) 
(G 58) 
(G 59) 
(G 60) 
(G 61) 
(G 62) 
!G 63) 
(G 64) 
(G 65) 
(G 66 ) 
(G 67) 
(G 68 ) 
(G 69) 
(G 70) 
!G 71) 
(G 72) 
(G 73) 
(G 74) 
(G 75) 
(G 76) 
(G 77) 
(G 78) 
(G 79) 
(G 80) 
(G 81 ) 
(G 82) 
( G 83) 
(G 84) 
(G 35) 
( G 86 ) 
(G 87) 
(G 88) 
(G 89) 
(G 90) 
(G 91) 
(G 92) 
(G 93) 
! G 94) 
(G 95) 
(G 96) 
(G 97) 
!G 98) 
(G 99) 
(G 100) 
(G 101) 
(G 102) 
(G 103) 
(G 104) 
(G 105) 
(G 106) 
(G 107) 
(G 108) 
(G 109) 
(G 110) 
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APPENDIX H 


SUBROUTINE VARLET 


THIS SUBROUTINE WAS WRITTEN IN CDC FORTRAN* VERSION 2.1* TO RUN ON CDC 6000 
SERIES COMPUTERS WITH THE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
MODIFICATIONS MAY BE REQUIRED PRIOR TO USE IN OTHER COMPUTERS. THIS SUBROUTINE 
HAS BEEN FOUND TO BE SATISFACTORY ON THE AFOREMENTIONED COMPUTERS. WHICH CARRY 
THE EQUIVALENT OF APPROXIMATELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRE- 
CISION MAY REQUIRE MODIFICATION TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS 
OF EQUAL ACCURACY. 


THIS SUBROUTINE IS REQUIRED BY THE PROGRAMS OF APPENDICES A, B, C, AND D. 



SUBROUTINE VARLET (ANGL) 

(H 1) 


COMMON ZETA, ETA , GAMMA , XOVERH , YOVERH , ZOVE RH * DELTA ( 28 ) 

(H 2) 


DIMENSION V( 1*9) ,ADEL(28) 

(H 3) 


SC = S I N( ANGL*0.01 745 329 25199) 

(H 4) 


CC = CO SIAN GL *0.01745229251 99) 

(H 5) 


Z6=ZETA*Z0VERH+1 . 

(H 6) 


Z8=-Z6 

( H 7) 


Z7=Z8-1. 

(H 8) 


DELTA! 5 ) =0.0 

(H 9) 


DO 20 M= 1 * 7 

(H 10) 


DO 20 N= 1 * 7 

(H 11) 


IF (N.EQ.4. AND.M.EQ.4) GO TO 20 

(H 12) 


DO 22 Jl = l * 9 

(H 13) 

22 

VC 1 * J 1 ) =0 . 

(H 14) 


ADEL ( 5 ) =0 . 

(H 15) 


AM=M-4 

(H 16) 


AN=N-4 

(H 17) 


X=ZETA*XOVERH 

(H 18) 


Y=ZETA*( YQVERH-2 •* AM*GAMMA+GAMMA* I l.-ETA)*! 1.— (-1 • )**M) ) 

(H 19) 


Z=ZETA*I ZOVERH—4 •* AN ) 

( H 20) 


A=SQRT(X*X+Y*Y+Z*Z) 

(H 21) 


B=A+Z*CC-X*SC 

(H 22) 


VI 1*1 ) = t I X*X+Y*Y )/( B*A*A*A) )-( IZ+A*CC )/< B*A ) )**2 

(H 23) 


Z=-Z-2. 

(H 24) 


A= SQRT ( X* X + Y*Y+Z *Z ) 

(H 25) 


B=A+Z*CC-X*SC 

(H 26) 


V( 1 * 3 > = ( ( X*X+Y*Y)/I B*A*A*A>)-( I Z+A*CC ) / ( B*A ) )**2 

(H 27) 


IF (ANGL. EQ. 90.0) GO TO 23 

IH 28) 


X=X-( SC/CC ) 

(H 29) 


Z=-Z-l. 

(H 30) 


A=SQRT(X*X*Y*Y+Z*Z) 

(H 31) 


B=A+Z*CC-X*SC 

(H 32) 


VI 1 *2 ) = ( (X*X+Y*Y )/{ 8*A*A*A) )-I IZ+A*CC)/( B*A) )**2 

(H 33) 


Z--Z 

! H 34) 


B = A+Z *CC— X *SC 

(H 35) 


V! 1 * 4 )-( ( X*X+Y*Y )/i 8*A*A*A) )-( I Z+A*CC ) / 1 B*A ) )**2 

{ H 36) 

23 

ADEL ( 5 )= I (-1. )**(M+N) )*( V ( 1, 1 )-VI 1 ,2)-V( 1 , 3 ) + V! 1 *4) ) 

(H 37) 


DELTA! 5)=OELTA(5)+ADEL (5) 

(H 38) 

20 

CONTINUE 

(H 39) 


X=ZE T A* XOVERH 

CH 40) 
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APPENDIX H 


Y=ZETA*YOVERH (H 41) 

Z=Z7 (H 42) 

A=SQRT(X*X+Y*Y+Z*Z) (H 43) 

B=A+Z*CC-X*SC (H 44) 

Vil,7) = ((X*XfY*Y)/( B*A*A*A ) ) - ( (Z+A*CC ) /( B*A) )**2 (H 45) 

IF ( ANGL.EQ.90.0 ) GO TO 26 <H 46) 

X=X-(SC/CC) (H 47) 

Z=Z6 (H 48) 

A=SQRT ( X*X+Y*Y+Z*Z ) (H 49) 

B=A+Z*CC-X*SC (H 50) 

Vt 1 ,6)=l ( X*X + Y*Y )/( B*A*A*A))-( (Z + A*CC )/< B*A) )**2 (H 51) 

Z=Z8 (H 52) 

B=A+Z*CC-X*SC (H 53) 

V( 1 1 8 ) = ( < X*X+Y*Y )/( B*A*A*A ) ) - ( ( Z + * *CC ) / ( B*A ) )**2 (H 54) 

26 DELTA! 5) = DELTA(5 ) — V ( 1 »6)— V( 1 1 7 ) +V ( 1 » 8 ) <H 55) 

AMT=-2.*GAM*A*ZETA*ZETA/3.141592E5358979 ( H 56) 

DELTA(5)=DELTA(5)*AMT (H 57) 

RETURN (H 58) 

END (H 59) 
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30 40 50 60 70 80 90 

Effective wake skew angle, X, deg 


(a) Effect of y. 0. 

Figure L- Effect of variations in width-height ratio and model height on vertical interference due to lift for a vanishingly small model 

laterally centered in a closed tunnel. 











Figure 4.- Sketch of high-speed section showing one scheme for folding ground plane to form tunnel walls. 




Effective wake skew angle,X,deg Effective wake skew angle, X,deg 


(a) Required width-height- ratio schedule 
for e w>L = 0. 


(b) Interference factors for variable and fixed 
width- height- ratio tunnels. 


Figure 5.- Required schedule and comparison of interference factors for unswept wings centered in a variable width- height- ratio closed-on- bottom-only tunnel, a = 0°. 
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Effective wake skew angle, X, deg Effective wake skew angle, X, deg 


(a) Required width-height- ratio schedule 
for 6 w,L = °* 


(b) Interference factors for variable and fixed 
width- height- ratio tunnels. 


Figure 6.- Required schedule and comparison of interference factors for unswept wings, laterally centered, but 0.1 H* below the centerline, in a variable width-height- 

ratio closed-on-bottom-only tunnel, a = 0°. 




Effective wake skew angle ,x, deg 


(a) Required width- height- ratio schedule for & w ^ = 0. 




30 40 50 60 70 80 90 

Effective wake skew angle, x, deg 


(b) Interference factors for variable and fixed width- height- ratio tunnels. 

Figure 7.- Required schedule and comparison of interference factors for unswept wings, laterally centered, but 0.2 H* below the centerline, 

in a variable width- height- ratio closed-on-bottom-only tunnel, a = 0° 
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30 40 50 60 70 80 90 

Effective wake skew angle, X, deg 


Figure 9.- Effect of model configuration on required schedule of width-height ratio for centered models in a variable width-height-ratio tunnel 

o = 0.5; a - 0°. 




Effective wake skew angle, X, deg 


Figure 10.- Effect of model configuration on required schedule of width-height ratio for laterally centered models mounted 0.1 H* below the 
centerline in variable width-height-ratio closed-on-bottom-only tunnel, o = 0.5; a = 0°. 



Effective wake skew angle, X, deg 


Figure 11.- Effect of model configuration on required schedule of width-height ratio for laterally centered models mounted 0.2 H* below the 
centerline in a variable width-height-ratio closed-on-bottom-only tunnel, a = 0.5; a = 0°. 
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Effective wake skew angle ,X, deg 


Effective wake skew angle, X,deg 


(a) Wing, A = 45°. 


(b) Rotor. 


Figure 12.- Effect of angle of attack on required schedule for centered models in a variable width-height-ratio closed-on-bottom-only tunnel, a = 0.5. 




Effective wake skew angle.x, deg 


(a) Wing, A = 45°. 



Effective wake skew angle, X, deg 


(b) Rotor. 

Figure 13.- Effect of angle of attack on required schedule for laterally centered models mounted 0.1 H # below the centerline in a variable width- 

height-ratio closed-on-bottom-only tunnel, o = 0.5. 
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Effective wake skew angle, X , deg 


(a) Wing, A = 45°. 



Effective wake skew angle, X, deg 


(b) Rotor. 

Figure 14.- Effect of angle of attack on required schedule for laterally centered models mounted 0.2 H* below the centerline in a variable width- 

height-ratio closed-on-bottom-only tunnel, o = 0.5. 
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40 50 60 70 80 90 

Effective wake skew angle, X, deg 


30 40 50 60 70 80 

Effective wake skew angle.X.deg 


(a) Interference factors at wing and tail. (b) Difference between interference factors at wing and tail. 

Figure 15.- Interference factors at tail in both variable and fixed width-height-ratio tunnels for centered models, o = o { = 0 ; a = 0°. 










Effective wake skew angle, x,deg 

(b) Horizontal interference due to drag. (c) Vertical interference due to drag. 

Figure 20.- Residual interference factors for a vanishingly small model centered in variable and fixed width-height-ratio tunnels. 
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Figure 23 .~ Model of variable model-height closed-on-bottom-only test section. 


L-67-10 000 






Effective wake skew a ngle,X J deg 

(a) Required schedule for 6 w ^ = 0. 



30 4 0 50 60 70 80 9 0 


Effective wake skew angle, X,deg 


(b) Interference factors for variable and fixed model height tunnels. 

Figure 24.- Required schedule and comparison of interference factors for unswept wings in a variable model-height closed-on-bottom-only tunnel 

having a width-height ratio of 2.0. 




50 60 70 80 9 1 

Effective wake skew angle, x, deg 


(a) Required schedule for 6 . = 0. 

W, L 



Effective wake skew aagle,X / deg 


(b) Interference factors for variable and fixed model-height tunnels. 

Figure 25.- Required schedule and comparison of interference factors for unswept wings in a variable model height closed-on-bottom-only 

tunnel having a width-height ratio of 1.5. 


Effective wake skew angle ,X,deg 

(b) Interference factors for variable and fixed model-height tunnels. 

Figure 26.- Required schedule and comparison of interference factors for unswept wings in a variable model height closed-on-bottom-only 

tunnel having a width-height ratio of 1.0. 





Effective wake skew angle, x, deg 


(a) Required schedule for 6 . = 0. 

IV, L 



Var £ COB 
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Effective wake skew angle, X, deg 


(b) Interference factors for variable and fixed model-height tunnels. 

Figure 27 .- Required schedule and comparison of interference factors for unswept wings in a variable model-height closed-on-bottom-only 

tunnel having a width-height ratio of 0.5. 






Effective wake skew angle ,X, deg 


Effective wake skew angle.X, deg 


<a) Interference factors at wing and tail. 


(b) Difference in interference at wing and tail. 


Figure 29.- 


Interference factors at tail in both variable and fixed model-height tunnels, y = 1.5; o - o t - 0; a - 0°. 




Effective wake skew angle ,x, deg 


Effective wake skew angle, X, deg 


(a) Interference factors at wing and tail. (b) Difference in interference at wing and tail. 

Figure 30.- Interference factors at tail in both variable and fixed model-height tunnels, y = 1.0; o = = 0 ; a = 0°. 
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30 40 50 60 70 80 90 30 40 50 60 70 80 90 30 40 50 60 70 80 90 

Effective wake- skew angle, x, deg Effective wake skew angle, x, deg Effective wake skew angle, x, deg 

la) Horizontal interference due to lift. (b) Horizontal interference due to drag. (c) Vertical interference due to drag. 


Figure 32.- Residual interference factors for a vanishingly small model in variable and fixed model-height tunnels having a width-height ratio of 2.0. 
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30 40 50 60 70 80 90 

Effective wake skew angle, X,deg 


(a) Horizontal interference due to lift. 



30 40 50 60 70 80 90 

Effective wake skew angle ,X, deg 


(b) Horizontal interference due to drag. 



30 40 50 60 70 80 90 

Effective wake skew angle,X,deg 


(c) Vertical interference due to drag. 

Figure 33.- Residual interference factors for a vanishingly small model in variable and fixed model-height tunnels having 

a width-height ratio of 1.5. 
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Figure 34.- Residual interference 
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Effective wake skew angle, X,deg Effective wake skew angle, X.deg Effective wake skew angle, X , deg 

(a) Horizontal interference due to lift. (b) Horizontal interference due to drag. (c) Vertical interference due to drag. 

Figure 35.- Residual interference factors for a vanishingly small model in variable and fixed model-height tunnels having a width-height ratio of 0.5. 
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(a) Var ? COB. 


(b) S = 1 closed. 


(c) Z = 1 COB. 


Figure 38.- Spanwise distribution of interference over an unswept wing in variable and fixed model-height closed-on-bottom-only tunnels having a width-height ratio 

of 1.5. a = 0.5; a = 0°. 



(a) Var £ COB. 


(b) £ = 1 closed. 


(c) ? = 1 COB. 


Figure 39.- Spanwise distribution of interference over a wing with 45° of sweep in variable and fixed mode I -height closed-on-bottom-only tunnels having a width-height 

ratio of 1.5. o = 0.5; a = 0°. 




(a) Var Z COB. 


(b) C = 1 closed. 


(c) ? = 1 COB. 


Figure 42.- Longitudinal distribution of interference over a rotor in variable and fixed model-height closed-on -bottom-on iy tunnels having a width-height ratio of 1 5 

o = 0.5; a = 0°. 







3 


r 

2.0 

, 5 

1.0 

.5 



0 I I I ... I I 1 1 1 1 1 1 1 1 

30 40 50 60 70 80 90 

Effective wake skew angle, X, deg 


(a) o = 0. 

Figure 43.- Relative recirculation of equal-area variable model-height closed-on-bottom-only tunnels. Reference tunnel has a width-height 

ratio of 2.0. Un swept wings; a = 0°. 
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30 40 50 60 70 80 90 


Effective wake skew angle, X, deg 

(a) Required schedule for A6 . = 0. 



Effective wake skew angle,X,deg 


(b) Interference factors. 

Figure 46.- Required schedule and resultant interference factors to reduce A6 W L to zero in a variable model-height closed-on-bottom-only 

tunnel having a width- height ratio of 2.0. a = = 0. 
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30 40 50 60 70 80 90 

Effective wake skew angle, X,deg 


(a) Required schedule for A6 W L = 0. 



30 40 50 60 70 80 90 

Effective wake skew angle,X,deg 

(b) Interference factors. 

Figure 47.- Required schedule and resultant interference factors to reduce A6 w L to zero in a variable model-height closed-on-bottom-only 

tunnel having a width-height ratio of 1.5. o = Oj = 0. 
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(a) Required schedule for A6 W L = 0. 
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(b) Interference factors. 

Figure 49.- Required schedule and resultant interference factors to reduce A6 w L to zero in a variable model-height closed-on-bottom-only 

tunnel having a width-height ratio of 0.5. a = <fc = 0. 









(c) Rotor, a = 0.5. 
Figure 50.- Concluded. 
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(a) Horizontal interference due to lift. 
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(b) Horizontal interference due to drag. 


(c) Vertical interference due to drag. 


Figure 52.- Residual interference factors for a vanishingly small model centered in a variable width-height-ratio closed tunnel. 
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(b) As a function of V/V*. 
Figure 57.- Concluded. 





in a variable model-height closed tunnel. 
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(b) Interference factors for variable and fixed £ tunnels. 

Figure 61.- Required schedule for minimum interference and a comparison of interference factors in a variable model-height closed tunnel. 

y = 0.5; o=0. 
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